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1 Introduction 



Both Spectral Geometry and Semi-classics seem to be concerned with more or 
less the same thing - studying the correspondence between classical-dynamical 
and geometrical quantities, such as the metric, the geodesic flow or the effec- 
tive potential for example, on the one hand, and on the other - the spectrum 
of a Laplace or Dirac-type operator, acting on the space of square-integrable 
functions on the given manifold. 

However it seems that each theory has its own independent motivation - the 
physicists' approach, referred to as 'semi-classical' is, widely speaking, a man- 
ifestation of the famous Bohr's correspondence principle^, while Spectral ge- 
ometry, representing the mathematician's point of view, is mostly dedicated 
to resolving inverse problems, or to the contrary, finding counterexamples. 
The essence of a typical inverse spectral problem is most naturally expressed 
by the question posed by Mark Kac in 1966, 'Can one hear the shape of a 
drum?', that is usually considered the birth of Spectral geometry. The idea 
is the following: given a smooth compact Riemannian manifold, one may 
construct, with the help of its metric, an unique Laplace- Beltrami operator, 
which is elliptic, self-adjoint, with discrete, definite spectrum, that obeys 
certain growth rate estimates^. It is then quite natural that the geometry of 
the underlying manifold determines the spectrum. The question is to what 
extend the same is true in the reverse order, or, in other words, can we re- 
trieve from the spectrum a sufficient amount of geometrical data in order to 
determine completely the metric or the shape of the boundary? It turns out 
that most typically the answer is negative. Even in dimension two there are 
examples of iso-spectral domains that are not isometric. Nevertheless, there 
are few cases of resolved inverse problems, such as low dimensional spheres 
and tori (d < 7), hyperbolic surfaces, real analytic surfaces of revolution and 
the result of Guillemin and Melrose, that is our main example. The methods 
for obtaining such results demand solid background in geometry, mechanics, 
PDE's and pseudo-differential calculus, as well as several more specific tools. 

1 postulating equivalence between the classical and quantum description in the limit 
h — > Q and certain analogy close to that limit - a good example of its validity is a famous 
theorem stating that the mean values of the quantum probabilities evolve according to 
classical laws 

2 in order to have all these properties we assume reasonable boundary conditions in the 
non-closed case and as for the dcfinitcncss, only a finite part of the spectrum can be with 
inverse sign in the Neumann case 
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That is why we dedicate the first half of the present work to the introduction 
of some preliminaries that are crucial for understanding the framework of the 
problem and apart from that give some physical insight. 
The second half already goes straight to the topic with the introduction of 
trace formulae, that appear to be a major tool both in Semi-classics and 
Spectral Geometry. These equalities relate infinite sums over closed classi- 
cal paths to traces of functions of the laplacian. Some of these are exact 
(Poisson and Selberg), while others (Gutzwiller, Berry-Tabor) rely on semi- 
classical asymptotic expansions. This common apparatus is one more reason 
to regard the two theories as being closely connected - one may treat the on 
the one hand semiclassical analysis as a tool of Spectral geometry, and on the 
other - consider a spectral-geometrical problem in the semiclassical context. 
The heat trace is somewhat exceptional - it expands the heat propagator 
as a power series in the high-temperature limit. As heat is diffusive, this 
expansion could not be related to classical hamiltonian dynamics, but is still 
rather sensitive to geometry - symbolic calculus yields infinitely many terms 
involving the metric derivatives when constructing a formal inverse and the 
coefficients appear to be linear combinations of the polynomial curvature 
invariants of the corresponding degree. This formalism appears to be very 
fruitful in rather unexpected directions. As we discuss in the following, it 
gives a regularization scheme and a handy tool for calculating conformal 
anomalies in quantum gravity, provides a simple procedure for obtaining the 
Korteweg-de Vries hierarchy and the infinite set of first integrals, but finds 
also an application in Index theory, statistical mechanics and others. 
Next we introduce quantum billiards as typical model systems in both the- 
ories. In spectral geometrical context they are considered as Euclidean do- 
mains in the context of the inverse spectral problem, while semi-classical (or 
geometrical-optical) approach treats them as quantum dynamical systems 
close to the classical limit. The most familiar interpretation is a particle 
captured in a box with infinite potential walls (Dirichlet boundary condi- 
tions). This setting, no matter how conceptually simple, appears to be a 
powerful tool for quantizing classically chaotic systems. Apart from this, 
spectral statistics can be used for retrieving qualitative information about 
the underlying dynamics, as we shall see in the following. 
In the last section we discuss some strategies and techniques for solving in- 
verse spectral problems and provide concrete examples in which they have 
found successful application. At the end we review the Guillemin- Melrose 
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problem^ and calculate the spectral variation with respect to varying the 
function in the homogeneous Robin boundary condition. The corrections, 
expressed as single-layer potentials, are found by using a slight modification 
of the standard in quantum mechanics apparatus of perturbation theory. 



3 The authors show that the function, participating in the homogeneous Robin bound- 
ary condition for the region in the inside of an ellipse is completely determined by the 
laplacian spectrum (more precisely, the resolvent trace built with its help) together with 
the spectrum of the pure Neumann problem 
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2 Dynamics and Geometry 



The classical dynamics of a physical system is usually derived from varia- 
tional principle. According to that formalism, we prescribe for example to 
a particle, constrained to move on a surface X and sensitive to external 
potential V, a Lagrangian function^ 

where (q, q) are the local coordinates in TX. 

With the help of C we easily construct the action functional 

S = / C{q, q) dr 

as s function over the space of paths with fixed endpoints go = q{to) and 
Qi — q(ti) respectively. The classical path is chosen as the critical point of S, 
that is the one for which the action is stationary, that is = 0. Substituting 
the above definition of S gives, after integration by parts 

Jt \~dq~) ~l)q~ U 

which is known as the Euler - Lagrange equation of motion. 
With our particular choice of Lagrangian, the Euler-Lagrange Equations cer- 
tainly coincide with the Newton second principle 

q = _ W(<7) 



2.1 Hamiltonian Dynamics 

Usually, however it is most convenient to describe such systems by means of 
Hamilton formalism. In order to switch from Lagrangian to Hamilton for- 
mulation of dynamics, we need to consider momenta p instead of velocities q, 
that is moving from the tangent to the co-tangent bundle TX — > T*X, q —> p. 
The Lagrangian C(q, q) undergoes a non-trivial change, known as the Legen- 
dre transform defined as 

H(q,p) =pq- C(q,q), q^p 

4 in general we define a Lagrangian density, but in the one-dimensional case, both 
constructions are certainly equivalent 
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The Lagrangian equations of motion, derived from variational principle are 
then equivalent to the canonical equations of Hamilton 

• - 9H ■ - dH (2) 

^ dq ^ dp 

where q and p are referred to as the generalized coordinate and momentum 
respectively. 

One of the reasons to prefer this representation is that it makes rather ex- 
plicit some nice properties of the conservative systems and is very often more 
convenient for various calculations, including quantization. This relief is due 
to the elegant framework of symplectic geometry we are allowed to use here. 
The key is that the cotangent bundle T*X, unlike the tangent, possesses a 
natural symplectic structure (a non-degenerate closed two-form), that in the 
above coordinates is given locally by u = dp A dq. Such simple form of to is 
called canonical and so is the pair p, q. The canonical pair is however not 
unique - we can always choose another pair of canonical variables in which 
the form uo looks in the same way. If the system is integrable, we can always 
find one pair J, 0, such that the hamiltonian depends only on the J's, but 
not on the 0's and thus the equations (J2]) look much simpler: 

j = = v (3) 

The first is obvious and in the second v is a constant, since the J - variables 
are all constants of motion and H depends only on them. 
Such variables are called action-angle variables and provide in many cases a 
powerful tool for integrating the hamiltonian system - each Jk reduces the 
phase space to the hyper-surface Jk = const and the final reduction gives the 
unique solution. Of course, one can choose J\ = H, since the 'total energy' 
of such systems is always preserved, and then follow a standard procedure 
for construction of canonical basis, analogous to the method of Gram and 
Schmidt in the metric case. 

If, however, one has found a full set of action variables J«, it is not evidence 
enough that the system is integrable: for example, they may appear to be 
functionally independent and the reduction in this case is not full. Another 
problem might appear if the reduction procedure is not commutative. There- 
fore we put the restrictions that the Jj's must be in involution with respect 
to the symplectic form, u*(Ji, Jk) = and that the differentials dJk must be 
all linearly independent. Under these conditions, the existence of the action- 
angle basis is already sufficient for a hamiltonian system to be integrable. 
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Sometimes it is more convenient to write the evolution equations using the 
Poisson bracket defined as follows 



WG} =^,^= £gg _*•££. ( 4) 

Using this formalism it is easy to express the equations of motion (j2J) in terms 
of this new structure as 

p={H,p} q = {H,q} (5) 

and in general, each physical quantity in a system with hamiltonian H is 
transformed with respect to the flow defined infinitesimally by x = {H, x}. 
Therefore we have that each quantity that commutes with the hamiltonian 
the system (with respect to the Poisson bracket) is preserved by the evolu- 
tion, defined by this hamiltonian. Such quantities are called conservation 
laws and the hamiltonian itself is one of these due to the fact that our struc- 
ture is skew-symmetric. Moreover, if we present another hamiltonian whose 
action preserves the initial one, then obviously the new hamiltonian is a con- 
servation law of the initial system (a version of Nother's theorem). 
Thus the problem of finding a solution of an integrable hamiltonian system 
is equivalent to the problem of finding a complete set of pairwise commut- 
ing linearly independent hamiltonian vector fields. In practice, however, it 
is often not an easy task and a lot of techniques have been developed to 
overcome technical difficulties. It is out of our aim to go into detail with this 
matter, but just as an illustration of a possible hidden trap for beginners: 
suppose that we know two of the conservation laws of an integrable hamil- 
tonian system of dimension n > 2 and we wish to generate the others. A 
straightforward method is to use the Poisson bracket - since this structure 
fulfils the Jacobi identity we have a series of commutators of the two pre- 
served quantities, that are preserved as well. Nothing, however, can be said 
about their functional independence, so there is some extra work to be done. 
In order to make things clearer, we may think of the dynamical equations 
(J2J) as of a vector field Xh, associated with the hamiltonian of the system 
and therefore, carrying the name hamiltonian vector field. We write instead 
of© 

X H {$) = u*(dH, d$) = {H, $} (6) 

The flow of Xh then defines the evolution of the hamiltonian system. In 
particular, this flow preserves the symplectic form as it is easy to be seen by 
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the expression for its Lie derivative^ 

L Xh uj = i XlI ° du + d o i XH u (7) 

The first argument is zero, since u is closed. The second one can be trans- 
formed to d 2 H = by using that on the flow, defined by H we have 

ix H uJ = dH 

that is one way to write the canonical equations. 

Since Lx H is a linear derivation, the hamiltonian flow preserves, together 
with u, all of its degrees up to n, which means all even dimensional pro- 
jections of the phase volume, including the volume itself [§. Therefore, it is 
not surprising that in hamiltonian systems asymptotic convergence is not 
possible (unlike Liapunov), nor are any types of attractors. All we can have 
are either integrable (or quasi-integrable) or chaotic, usually ergodic types of 
motion. 

Another footnote to make here is that the notion of a hamiltonian vector 
field makes the analogy between Poisson brackets and commutators far more 
explicit. Namely, if a given quantity, say F, commutes with the hamiltonian 
H of the system (in the sense that {F, H} = 0), it may be regarded as a 
hamiltonian associated to a flow 

F — > Xp — > (p s 

that leaves our initial system undisturbed, which can be written as [Xjj, Xp] = 
and vise versa: each time the latter condition holds, we have a preserved 
quantity F for Xjj (or H for Xp) which can be found by exploiting the 
differential of the action, known as the integral invariant of Poincare-Cartan 

dS = pdq- Hdt (8) 

Since, on the one hand, <p s is a flow that leaves H and therefore its Legendre 
transform - the Lagrangian invariant, and on the other t— and s— derivatives 
commute (due to the vanishing of the above commutator), we may write 

q d dS ^ { ' H) d f ,dq jjdt\ 
ds dt ds dt \ ds ds J 

5 this formula, due to Cartan gives the derivative of a quantity in direction to the flow, 
defined by Xh (Lie derivative of a vector field). Here ix H is an operator, fixing the first 
argument equal to Xh 

6 note that u n ~ dVol(X) since dirndl" (X) = and note also that this imposes one 
nontrivial restriction on the field Xh, namely divXn = 
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and the expression in brackets on the righthand side is a constant of motion. 
Let us go back to the symplectic structure u>. Note that it splits the 2n- 
dimensional phase space into two subspaces on which it vanishes (in our 
initial basis it were the coordinate base and the cotangent momentum sub- 
space). When we change coordinates the splitting may not be so simple, 
but each time we have an n-dimensional submanifold, the restriction of uj on 
which is identically zero (take for example the submanifold, defined by the 
constants of motion J&) we call it a Lagrangian submanifold. 
Lagrangian submanifods, as we shall see shortly, are crucial both for dynam- 
ics and geometric optics. One way to construct such manifold is to start with 
one that we know (say the g-space) and make it flow according to a canonical 
change of variables. The latter means simply that u should be preserved by 
the so chosen flow. The group, that is responsible for such canonical trans- 
formations is called real symplectic group and denoted usually by Sp(2n, R). 



Generating functionals and Hamilton-Jacobi equations: In practice, 
the canonical change of variables, we are interested in, is being made via gen- 
erating functional. The main idea roots in the observation that a coordinate 
transition <p : (p, q) —>■ (P, Q) is canonical if the co-homological class of the 
momentum one-form is preserved under this transformation, that is 

pdq = PdQ + d& 

One can easily assure oneself by differentiating both sides of the above equal- 
ity that, since gP$ = 0, is volume - preserving and hence, canonical. The 
potential $, expressed in (q, Q)-variables3, is called a generating functional 
for the transformation and appears to be a useful tool for the description 
of the canonical equations in the new variables. Moreover, the flow induces 
its pull-back, acting on the hamiltonian functions <fi* : H — > H which can 
be used as an alternative description of the theory. In particular, we may 
use this pull-back to find the generating functional of a transition to action- 
angle variables and then express the hamiltonian in these new variables - that 
would simply be H(Q). Since, by the definition of a generating functional 

7 this can always be done as long as the Jacobian det — — is non- vanishing, or, equiva- 

op 

lently, the Hessian of $ is non-degenerate 
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we have 

_ <9$ _ <9$ 

our aim is to find a solution of the non-linear PDE 

H(Q) = H(—,q,t) (10) 
known Hamilton- Jacobi type equation. 

The action variables Q are the constants of motion defined by p = -^—(Q, q) 

dq 

and the linear flow on the 'angles' is naturally given in terms of the new 
hamiltonian: 

P = JJL 

dQ 

Thus, resolving the above non-linear PDE is equivalent to finding the tra- 
jectories of the initial Hamiltonian system. This equivalence can be utilized 
in both directions, but, surprisingly enough, it turns out to be the most pow- 
erful tool known so far for integrating the canonical Hamilton equations of 
motion. 

At the end we mention that the classical action also satisfies a Hamilton- 
Jacobi type equation. Namely, it is straightforward to check that 

d t S + H( — ,q,t) = (11) 

due to the explicit expression for the Poicare - Cartan integral invariant (jSJ). 
As we shall see shortly, this leads to an analogy between the description of 
the dynamics of a classical particle and that of phenomena from geometric 
optics. 



2.2 The Geodesic Flow and Geometric Optics 

A good example of a hamiltonian system can be seen in the problem of a 
free particle, constrained to move on a Riemannian surface M. It is well 
known that all possible trajectories are namely the geodesies on M. Thus 
the concept of a geodesic flow is somewhat central for classical mechanics, 
let alone that even if we have interaction in the system, its potential might 
be 'transformed' into some effective curvature of the underlying manifold, so 
at the end we have again the problem of a free particle. This is the case 
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with the two-bodies problem in physics - the solutions are easy to obtain as 
geodesies over the sphere. 

To put these thoughts in order, let us first consider the metric Lagrangian 
L = - (q, q) , which is nothing but the lagrangian of a free particle of unit 
mass, constrained to our manifold. The corresponding hamiltonian is then 
given by H = -(p,p). This immediately yields conservation of momentum 
and extremal trajectories. 

Another way to find the equations of geodesies is using the concept of co- 
variant differentiation on M which is locally given by 

(W), = v + V) h v k 

where the additional term accounts for the transformation of the local basis: 
Y l - k = dje\ is called the Lie-Christophel Symbol. With a suitable choice of 
the evolution parameter the equation of the geodesies looks like 

v i + r i jk v j v k = (12) 

Both formalisms are certainly equivalent and we have the choice between 
the former Lagrangian (or Hamiltonian) representation, fully in the spirit 
of extremal principle and the latter, which accounts for the fact that the 
tangents to the geodesies are auto-paralle@. 

Let us consider now a wave function with rapidly varying^ phase <3>. When 
substituted into the Schrodinger equation, semi-classical approximation gives 
up to first order in h a Hamilton - Jacobi equation in the form 

d t $ + H(— ,q,t) = Q (13) 

This means, on the one hand, that semiclassical 'trajectories' of photons are 
given by Hamiltonian canonical equations with principal function $ and, 
on the other, these trajectories are all perpendicular to the wavefront $ = 
const (Huigens principle). There is therefore a kind of dualism between the 
description of geometric-optical phenomena by means of their wavefronts, or 
by the light rays. 

In the present considerations, we pay much more interested to the second 



trivially translated along their own integral curves 
9 compared to the alterations of the potential, or the geometry 
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interpretation, as leading to more obvious optical - mechanical analogy. For 
example, by Fermat principle, light trajectories minimize the time of arrival, 
which is proportional to the functional known as optical distance 



where n is the refraction coefficient of the media and dl - the length element. 
In particular, for homogenous media (n = 1), those trajectories are nothing 
but the geodesies. For example, it is well known from general relativity, 
that on the large scale, the light travels through the Universe, following the 
space-time geodesies with respect to the metric, 'curved' by the presence of 
matter, according to the famous Einstein - Hilbert equations^. 

2.3 Morse Theory and Catastrophes 

We take as our main example the geodesic flow or, alternatively, the geomet- 
ric optical problem concerning ray families. Both descriptions happen to give 
somewhat equivalent picture due to Fermat principle of stationary action. 
To begin with, let us recall the basics from the classical Morse theory. The 
first crucial statement of the theory is that on each compact Riemannian 
manifold M there exists a smooth function / with only simple, and thus 
isolated, critical points^], we refer to as Morse function. Then at each such 
critical point the / is locally isomorphic to a quadratic form, so that local 
extrema are given as minima, maxima and saddle points only. Define the 
index of this quadratic form, called the Hessian form, as the number of its 
negative eigenvalues. Then the main theorem of the classical Morse theory 
states that if to each critical point Xj we assign a cell of dimension equal 
to the index of Hess(/) at Xj and then 'glue' all the cells starting from the 
lowest dimensional we obtain a smooth manifold with the same homotopy 
type as M. 

10 stating that the Riemann curvature tensor of space-time is proportional to the stress- 
energy tensor of the classical fields (which is, on the large scale, mostly dominated by 
gravitation) 

n the point xq is called critical for V/ = at xq and simple, or non-degenerate if in 
addition to that the Hessian VV/ is non-degenerate at this point 

12 the boundary of each cell is glued to the cells of lower, but maximal possible dimension: 
for example a 1-cell glued to a 0-cell is isomorphic to a circle, a 2-to-l-gluing gives a disk 
and 2-to-0 makes a 2-sphere etc. 



S = c 
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The simplest examples usual, the two-dimensional sphere and torus, 

with the height function. For the sphere we obviously encounter a 0-cell 
(for the minimum) and a 2-cell (for the maximum) which is all we may ever 
expect. On the torus (if we make it 'stand' properly) we have two saddles 
plus the above ones, which gives two circles attached to a common point and 
then 'dressed' in a 2-cell skin. 



Back to the geodesic flow: Now consider the manifold of all piecewise 
smooth curves over a Riemannian manifold M and denote it by Q(M). 

Now functionals over M are recognized as functions in Q. In particular the 

1 f T1 

geodesic action S = - ||7(r)|| 2 (ir is one such function and we are going 

2 J T 

to pay special attention to it in the rest of this section. 
The variation of the classical functional £(7) is generated by vectors in T 7 f2. 
One such vector is the velocity field v — 7 along the path 7, or we may choose 
a transversal direction, say u G T^Q. The critical points of the functional 
with respect to one such variation give the stationary paths of a free particle, 
bound to move on M. The Euler-Lagrange equation defines the geodesic con- 
necting 7(tq) with 7(71). This geodesic however does not have to be unique 
- what happens in practice is that uniqueness is fulfilled naturally for close 
enough initial and final point, but globally it is no longer so. For instance 
there are infinitely many geodesies, connecting the south and the north pole 
of a round sphere. In order to understand this phenomenon one has to fix a 
geodesic and take into account the second variation of S around it as well. 
Uniqueness is required whenever the Hessian^! 

5lS(l) = ~ r (v, + i?( 7 , u)^)dr (14) 

J TQ 

is strictly definite, but as long as one of the eigenvalues changes sign this 
uniqueness is broken. To be just a little bit more precise, at the moment 
of the sign transition the Hessian is degenerate, meaning that the field u 
satisfies the Jacobi equation 

V^x + R(j, u)7 = (15) 

13 for derivation of this formula one may refer to the famous Milnor's comprehensive 
book 'Morse Theory' 
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The field u is called a Jacobi field and it deforms the initial geodesic 7 into a 

smooth one-parameter family of geodesies g(s, r). We just outline the proof 

which comes together with the proof of the above formula f|T5|) . 

The first variation of the action gives the Euler - Lagrange equation of the 

form 



since the lagrangian is independent of coordinates^. Then, varying second 
time in the w-direction (u, v G TQ(M)) and assuming that both variations 
are independent, we simply arrive at the second order equality 



provided [u, v) = 0, which can be easily arranged at least locally Secondly, 
due to equality of the second derivatives, we have V„7 = V^u and thus we 
easily recover ffTBT) and besides, if one seeks for a smooth deformation leav- 
ing the geodesic property of 7 invariant, that is V^g(s, r) = preserved in 
the M-direction, after a while one comes upon the requirement that u is a 
Jacobi field and vice versa - to each Jacobi field at certain geodesic 7, one 
may associate a one-parameter smooth deformations, preserving its geodesic 
property. 

Now, given a geodesic 7, two points on it are called conjugate if there exists 
a Jacobi field on 7 which vanishes at both points. The number of such lin- 
early independent fields is called multiplicity of this conjugacy. This number 
is always finite since it must be less than the dimension of the underlying 
manifold. Then the Morse index of the geodesic 7 is defined as the number 
of conjugate points along 7, counted with their multiplicities and this num- 
ber is equal to the index of the quadratic form VVS" at the critical point 
7. Furthermore, the main result of the so extended theory states that chord 
space of geodesies with given endpoints is homotopy equivalent to a union 
of cells of dimensions, equal to the Morse indices of all possible geodesies 
having the same endpoints fixed. 

14 this is easy to derive by fixing the endpoints and integrating by parts 





Now, we have to keep in mind that 



V u V^ 7 (r) = V^V U + R(u, 7) 
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First we note that the way we defined the Morse index of the geodesic flow, 
it happens to coincide with the Maslov index of the Hessian of S, defined 
as the number of sign changes of the eigenvalues from negative to positive 
minus the number of opposite changes. 

Second, we wish to discuss in brief the Jacobi equation (TT51) . It is a second- 
order evolution equation that describes the behavior of the geodesic flow on 
M. Let us take the two-dimensional compact surfaces as an example. The 
constant scalar curvature allows for writing (fl5l) in the form u + ku = 0, 
where k is the scalar curvature and u is the vector of deformation chosen 
to point in a direction normal to the one defined by the evolution param- 
eter. Then for the sphere the solutions are oscillating, which is a hint for 
the presence of conjugate points^. On the torus, as it should be expected, 
the geodesies are parallel, so the flow is uniform, while on hyperbolic sur- 
faces it is exponentially divergent. This is a nice illustration of the fact that 
whenever the sectional curvature in all two-directions is strictly negative, all 
Morse indices are zerc0. 

The theory has two quite natural extensions - one of them is to consider 
closed orbits, but them one needs to be careful about periodic conditions 
imposed on the deformations and possible complications due to revolutions 
of simple periodic orbits. The results, however, are quite analogous to what 
we have seen so far. 

Another generalization is to consider manifolds with billiard types of geodesies, 
where j(r) may suffer noticeable jumps. In this case there are additional 
terms in the expressions for the action variations, proportional to these jumps 
(and their square). Since the results are quite analogous and nothing con- 
ceptually new is obtained, we omit these considerations here, although we 
shall need the notion of Morse index namely in this case. 
Although a lot of progress has been done so far in Morse theory - there are 
generalizations to gauge and super-symmetric theories, we modestly refer to 
the basics, given in [39]. For Hamiltonian mechanics and symplectic geome- 
try, a good choice would be [1] together with the appendix. The Hamilton- 
Jacobi method can also be found there. 

15 more generally on S n the two poles are conjugated with the maximal possible multi- 
plicity n— 1 

16 this result is far not difficult to prove together with its almost immediate consequence 
that closed geodesies on closed manifolds of negative curvature are isolated, so that there 
is only one in each homotopy class - a simple way to do this relies on Bochner technique 
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The Poincare map Now is the time to briefly introduce the Poincare 
map for periodic orbits. Given one such orbit 7 with period T, so that 
7(T) = G T 7(0) = 7(0), one takes a point 1G7 and considers a symplectic 
transversal to 7 at x. Then the Poincare map is defined as the first return 
map to S 1 : 

P.-.s^s,, P 7 (C) = G T ^(C) 

where T(() is the smallest time it takes ( to return to S 1 
The linear Poincare map is then P 7 = dP 1 . Sometimes it is convenient 
to define P 7 on the span of complex normal Jacobi fields along 7. It is a 
symplectic vector space with respect to the Wronskian 

w(X,Y) = g(X,D a Y)-g(D a X,Y) 

Then P 7 is defined on this space simply by 

P 7 F (r) = Y(t + T) 

Note that the abstract definition we gave for P 7 is often reduced to rather in- 
tuitive constructions. For example, in the case of convex billiards, we take the 
symplectic transversal (or Poincare section&s they usually call it) to be the 
domain boundary, fibred with the normalized momentum (in more math- 
ematical terms the unit ball bundle BdQ) and then one defines P 7 as the 
first-return map of a neighborhood of an initial point £0 in BdVL. It could be 
interpreted also as an action of monodromy around a given closed geodesic. 
More precisely P 7 maps the variation 5£ = (8%, 8po) into 5£ relative to a 
circuit around the closed path 7. 

One major advantage to describe the billiard flow by means of the Poincare 
map is that in this way one reduces the n-dimensional dynamical problem 
to the study of a discrete iterative map on the boundary. As we shall see in 
the following, this is a great relief. 

The eigenvalues of the linear symplectic map P 7 help us to distinguish be- 
tween different types of closed geodesies. For example, 7 is said to be non- 
degenerate, if det(id — P 7 ) 7^ or, in other words, unity does not belong 
to Spec(P 7 ). Such non-degenerate closed orbits fall in one of the following 
classes: 

1. elliptic if all eigenvalues are of modulus one, so they come in families 
e ±ia \ ctj eRj'=l..n 
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Figure 1: The evolute of the ellipse is built as the envelope of all normals and is 
thus a caustic for the normal fibration 

2. hyperbolic if all eigenvalues are of the type e ±Mj , /ij G K. 

3. loxodromic if Xj = e ±ia ^ ± ^ 

Since P 7 is symplectic, no other possibilities are available. 

Lagrangian singularities So far we managed somehow to avoid the ques- 
tion of degenerate singularities both in our geometrical and dynamical con- 
siderations. This is the point where Morse theory fails and we need new tools 
to investigate in depth such structures. These structures themselves carry 
the curious name catastrophes and have been somewhat unrevealed mystery 
of geometry until the second half of 20 th century when mathematicians at 
last came up with a consistent theory on the matter. 

To illustrate better the idea of what a catastrophe might be, let us first give 
an example. Consider an ellipse in the plane and a nice Morse function on it, 
which will be the distance to a fixed point in the inside region. Now suppose 
that, by some unlucky coincidence, one has chosen a point which is a center 
of curvature. In this case the distance is no longer Morse function, since the 
ellipse looks like a circle up to second order from one such point of view. As 
it is well known, all centers of curvature of an ellipse form a piecewise smooth 
curve, which is also the envelope of all lines normal to the boundary, that is 
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Figure 2: The simplest catastrophe - the fold in two dimension, obtained by 
projection 

called the evolute of the ellipse. 

As we may see from this example, unlike Morse critical points, degenerate sin- 
gularities need no longer be isolated, but could come in various co-dimensions 
instead. 

Another typical example of a catastrophe are the caustics in geometric 
optics. Caustics are, widely speaking, envelopes of ray families and they 
are usually visible as a projection of focused beams. The first example that 
comes to mind is the focused through the reading-glass sunbeam that burns 
the paper. But this type of singularity is highly unstable - we can easily 
distort it by a small perturbation. Then it starts to look like a cusp - the 
curve from the previous example. 

All classes of topologically stable caustics, however, at least in low dimen- 
sional, have been studied and can be locally represented by polynomials. 
Let us become familiar just with the first few of them - the rest can be easily 
found in almost every book or article on catastrophe theory. 
In the maximal codimension we have only two possibilities - the fold, locally 
given by the potential function x 3 + ax. The parameter a is responsible for 
the 'unfolding' of the singularity. 

17 more precisely, dimensions up to twenty-two and co-dimensions - up to seven 
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Figure 3: The cusp catastrophe is the second and last possibility for structurally 
stable Lagrangian singularity in dimension two (codimension one) 

Common examples of a fold catastrophe are the sea's sparkles or the ap- 
pearance of rainbow via complete reverse refraction of light in small drops of 
water. 

A little more complicated is the cusp, which looks like x 4 + ax 2 + f3x. We 
already have two good examples of a cusp catastrophe, and the picture below 
shows how this singularity is projected in the plain. 

The swallowtail has co-dimension two and the corresponding polynomial is 
x 5 + ax 3 + f3x 2 + jx. For examples of the swallowtail and the more compli- 
cated catastrophes in practical problems we refer to [66J. 
As an easy exercise, one may check how variations of the control parameters 
alter the extrema of the potential function. There are critical values that cor- 
respond to bifrucation of different types: collision and separation of minima 
and maxima (and in particular anihilation of a minimum and a maximum), 
loosing stability (minimum to saddle point) etc. 

The whole catastrophe theory classification we talk about is somewhat spe- 
cial, although it includes almost all possible topics of interest both in ge- 
ometry and physics. Its peculiarity originates in the fact that we are con- 
sidering Lagrangian catastrophes only. We remind that a submanifold of a 
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2n-dimensional manifold is called Lagrangian if its dimension equals n and 
the symplectic form vanishes on it. The most trivial case is the base of a 
cotangent bundle, considered as a symplectic manifold (in particular - the q- 
plane) or, as it was in our geometric optical settings, the level set of the phase 
function or the integral surface of the Hamilton - Jacobi equation. Then we 
naturally call Lagrangian the singularities between lagrangian manifolds of 
equal dimension. In our case these were the projections of critical points on 
the coordinate space. This locus of branch points plays the role of an in- 
variant curve (or surface, or point) to which infinitely many rays (or particle 
trajectories, as far as classical dynamics is involved) tangate. In the case of 
the elliptic billiard table for example, the caustics are namely the confocal 
ellipses and hyperbolae with the help of which one integrates the system. 
An alternative treatment would pay interest in the wave-front singularities, 
instead of those of the light rays, but we are not going to consider this case 
here. The reader may look in [3J and [12] instead. In the appendices of [TJ 
one may find a lot of interesting considerations concerning geometric optics, 
normal forms and catastrophes. 
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3 Elementary Quantum Mechanics 



In section one we described the dynamics of a particle by a hamiltonian 
function together with the symplectic structure (or Poisson bracket) in the 
cotangent bundle. Then integration in this formalism, reduces to finding a 
complete set of commuting first integrals. 

In elementary quantum mechanics, as we shall see in an instant, this picture 
is somewhat preserved: the reduction of a stationary quantum mechanical 
system for example, which is this time represented by a linear operator, act- 
ing in Hilbert space, is equivalent to finding a generator of symmetry for the 
system. Since the corresponding operator commutes with the hamiltonian of 
the system, it provides a convenient 'splitting' of eigenstates. Various prob- 
lems have been attacked in this way, among which the hydrogen atom, the 
oscillator, the Zeeman effect and many others. A good comprehensive book 
on the basics is [23]. However, for the path integral formulation we refer to 
[27| and PB]. 



3.1 Heisenberg and Schrodinger Pictures 

We start with a picture in which each physical observable Z is represented 
by a linear self-adjoint operator, acting in the Hilbert space L 2 of square 
integrable functions. The evolution of a physical observable given again in 
terms of a hamiltonian and a Poisson bracket 



but this time the hamiltonian function is an operator itself and the bracket is 

given by a commutator times a numerical factor — with h being the famous 

in 

Plank constant. 

Straightforward integration of the above formula gives for the flow of Z 



The above expression is an adjoint action and expands in an infinite series 
of commutators according to the famous Campbell - Housdorff formula. 

IT 

It is sensible to call the operator U(r) = exp [— — H] the propagator associ- 
ated to H. Since H is self-adjoint by definition, the propagator is always 



Z = {H,Z} 



(16) 




(17) 
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unitary. Of course, H may depend on time itself and one should consider 
T-ordered exponents instead, but for the time being we neglect this case. 
The reason we have to call an infinite dimensional operator a physical observ- 
able is that what we can really observe and measure in experiments are the 
eigenvalues of the operator. Therefore we have to choose a canonical basis 
in which the operator is diagonal. One may think of the process of observa- 
tion or the experiment itself as being equivalent to one particular choice of 
such basis. But we know from simple algebra that two linear operators share 
a common canonical basis if and only if they commute. This is the place 
we come across the first crucial difference between classical and quantum 
mechanics, known as the Heisenberg uncertainty principle: in the quantum 
picture we are 'allowed to observe' only quantities that commute with each 
other. In deed, as we shall se in the following we are not even able to mea- 
sure simultaneously the initial coordinate x and the momentum p, so what 
is usually referred to as quantum indeterminism is laid in the very roots of 
the theory. 

Nevertheless, the equations of quantum mechanics give rather exact and in 
many cases the only known description for various phenomena from the world 
of particles, atoms and molecules, to the description of large ensembles such 
as polymers, crystals, gases and even stars. In order to obtain these equa- 
tions, however, we need to start with choosing a suitable basis. 
Let {4>k} be the basis of eigenfunctions for the observable W. We have 

W(j) k = w k 4>k 

and since W is self-adjoint, all eigenvalues are real and the 0j's are orthogonal 
with respect to the standard scalar product in L 2 : 



as long as they do not fall in the same degeneracy subspace (have the same 
w as an eigenfunction). 

We use the above to express the matrix elements of another observable Z as 

z ki = | Z4>i) 

Moreover, let ip = ^c^k be an arbitrarj@ unit vector in L 2 . For reasons 
that were already explained, we call the vectors in L 2 states. Then, in the 

18 once more, for the sake of simplicity we avoid the question of continual indices, but 



there might be an additional term of the form / c(k)<fi(k)dk in the generic case 
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state ipi the observable Z projects as 



(Z)^ = (ij\Zij) 



In view of the above definition of the scalar product and the normalization 
of ip, (Z)\$ can be regarded as an expression for the expectation value of Z 
in the state ip. With this in mind it is easy to interpret \ip\ 2 = as a 
probability density and the Fourier coefficients in the above expansion give 
respectively the transition amplitudes \ck\ 2 - the probabilities for transition 
from a state ip to fa, or the probability to detect Wk in the experiment. 
The above description (known as Heisenberg picture or matrix mechanics) 
gives a good intuitive link between Hamiltonian dynamics and quantum me- 
chanics. For practical reasons, however, it is often convenient to consider a 
dual situation, in which the observables remain unchanged in time, but states 
evolve instead. In explicit terms it follows directly from the above that if we 
assume Z{r) = Z(0) = const and ip(r) = U(r)ip(0), we easily obtain the 
same values for the 'measurable' quantities - eigenvalues, expectation values, 
transition amplitudes etc. The latter leads immediately to 



known as the (non-stationary) Schrodinger equation. Now it is time to reveal 
the explicit form of the operator H. In fact, just as in the classical case, it 
is a sum of kinetic and potential term, usually in the form 



agrees with the classical one, since obviously 

[p, x] = ih 

Now it becomes clear why coordinates and momenta are not compatible in 



19 in order to be fair, we add that this choice is only a convention - one may alternatively 
choose p — > p. and x — > —itV\7 p - the so called 'momentum representation' 
20 note that p and x are related via Fourier transform 




(18) 
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Thus the more precise we measure the coordinate, the less we know of the 
momentum and vise versa. The same remark certainly refers not only to p 
and x, but to any other pair of operators that satisfy analogous commutation 
relations, e.g. energy and time, angular momentum and angle etc. 
Let us now return to the equation ffT8]) . It is a second-order differential equa- 
tion that in most cases is practically unsolvable. Moreover, regarded as an 
evolution equation, the Cauchy problem for which does not exist (due to the 
uncertainty principle just discussed), it resembles more the wave equation, 
than the equation of conductivity, thanks to the i-factor in front of the time 
derivative. That is why its solution is sometimes called wave function. In 
many cases the solutions really behave like waves. 

Now consider a special class of solutions, for which the probability density 
does not evolve in time. This requirement allows ip to depend on r only 
through a factor of unit modulus, that is 



The eigenvalues of the Hamiltonian operator represent the energy spectrum 
Ek = fiuJk of the system. The spectrum is discrete whenever there is an 
appropriate energy bound. In the generic case we have an discrete and a 
continuous part (for example in scattering problems). 

The representation in which the Hamiltonian is diagonal is called therefore 
the energy representation and the above equation that defines it - the sta- 
tionary Schrodinger equation. We are going to use it a lot in the following. 



Let us consider one of the simplest classical problems and try to quantize. 
The classical hamiltonian of the linear harmonic oscillator (with unit mass) 
is given by 



where u is a constant of elasticity. If we apply directly the formalism we have 
built up till now, we the quantum hamiltonian in the form of a second-order 




to a Helmholtz type equation for the function i/j(x): 

= hui> (21) 



3.2 The Oscillator 
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differential operator with a quadratic potential term. We are interested in 
the eigenvalue problem for this operator. 

Following an idea of Dirac, we express the hamiltonian in terms of the oper- 
ators 



(Q + iP), a! = -=(Q-iP) (22) 



<2uj V2uj 
where P and Q are the quantum momentum and coordinate operators, re- 
spectively. 

It is easy to see that a and a) are not self-adjoint, but pair- wise adjoint, 
instead. Moreover, they fulfill the commutation relation {a, a^} = 1, that 
gives, by induction {a, {a)) k } = k(a^) k ~ 1 . In this way a acts on the functions 
of a) as formal differentiation. Now the quantized hamiltonian can be written 
in the form 

H = uo ( aa) I — u [ aa) H — 

V 2) \ 2 

which allows for computing easily 

[if, a] = —uja, [if, a)} = uoa) 

This construction already resembles very much the standard weight decom- 
position of semi-simple Lie algebras. The analogy goes a bit further with 
finding the highest (lowest) weight. 
More precisely, the spectrum of H is bound from below. 

E > — 
~ 2 

To see this, one simply needs to exploit the standard scalar product in L 2 
for an eigen-state, corresponding to the energy E. From the above, we have 
directly 

(^ E \H^ E ) = E\\i> E \\ 2 = uj\\a^ E \\ 2 + |||^|| 2 

which proves our assertion and gives in addition that equality is reached if 
and only if aip = 0. As we shall see shortly, such if) exists and E = — is the 

lowest state energy, also referred to as the 'vacuum energy'. 
In order to show that this lowest energy is reached, we consider the series 
of vectors ipk E -m = a m ipE- We claim that if ipE is an eigenvector with 
eigenvalue E then ip kE _ m is also such with eigenvalue E — mu. For m = 1 
we obviously have 

HaipE = aHipE — ojipE — (E — u)ipE 
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and the rest follows by induction. Since this series is diminishing and bounded 
from below, it should be finite, so there is such a natural number /c^ for which 
a kE+l ipE = and correspondingly such tp for which 

aipo = 0, Hi/} = -ip 

On the other hand, the operator a) generates infinite, non-vanishing series of 
solution tp n = (atifto) with eigenvalues E n = uj + respectively. 

Now it becomes clear why a* and a are called creation and annihilation 
operators respectively. With their help we can construct explicitly the basis 

of solutions. We start with ip and substitute in the differential equation 

00 

E — — . We obtain 

where the constant C is defined by the normalization condition. 

For the rest of the eigen-states, we obtain, by a complicated and non-interesting 

procedure 

ij n ~ H n (y/iux)e~^~ 
where H n are the Hermit polynomials. 



3.3 The Density Matrix 

Sometimes, when we are interested in the practical side of the matter, espe- 
cially in the study of ensembles in quantum statistics, it is more sensible to 
study directly partition functions (probabilistic distributions) for a certain 
problem, rather than the 'non-physical' wave-function ip. As we already saw, 
the meaningful quantity that we may get out of ip is its squared modulus 
ip*ip. What we do in practice, however, is to define instead the projector 
over the ip - state IL/, as acting on a state a by 

li^a = (ip\o-)ip 

and in Dirac notations we write n^, = \ifj) Note that it is in fact an 
integral operator with unit trace since, by normalization, = 1. This 

very projector we call a density matrix for the state ip and denote by p or 
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The reason for choosing this name is somewhat obvious: once a complete 
orthonormal basis of pure states fa is introduced in our Hilbert space the 
above integral operator adopts indices, namely 

Pkl = c* k c h Cj = {<f>j\ip) (23) 

and the diagonal elements Wi := pa can be regarded as probabilistic weigths 
of the 'pure states' fa. If the system is entirely in a pure state, then only one 
of them equals unity (all other matrix entries vanish) and p 2 = p. 
With the help of the density matrix it is easier to express statistical data in 
covariant terms. For example, the expectation value of an observable Z is 
now written as 

(Z) = tr(pZ) (24) 

and various thermodynamical functions can also be given in this way. For 
example the free energy F = —kTlntrp, or the entropy S = — ^ u^lniOfc. 
Note that this picture is still Schrodinger rather than Heissenberg and due 
to equation (fT8l) we have for the evolution equation for the density matrix in 
the form 

p = {p,H} (25) 

which resembles the dynamical equations for observables in Heissenberg rep- 
resentation, except for the 'minus' sign. 

In statistical mechanics, we know by general considerations that, up to a 
normalization factor, the diagonal elements of p are given in the energy rep- 
resentation by 

w k = e~P Ek 

where (3 is proportional to the inverse temperature^] and Ek is the energy of 
the corresponding pure state. In a more covariant form, we may write 

p = exp \—(3H] 

which leads directly to the famous Bloch equation for the density matrix 

-^P + Hp = (26) 

This is a diffusion or heat equation as it is easily seen from the simplest case 
- a free particle in Euclidean space. Then the solution is easily seen to be a 
Gaussian distribution. 



21 equals the inverse absolute temperature times the inverse Boltzman constant, to be 
more precise 
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3.4 Path Integrals 

Now we shall develop an alternative approach to the study of the density ma- 
trix. The so-called path integrals involve the quantum mechanical concept of 
non-locality, by letting the present state of the system be influenced not only 
by a single evolution, but instead, by all possible paths that link the initial 
and final points in phase spaced. Following Feynman, we write formally 



p(x in ,x f ,u) = exp 



Hu 



which holds for sufficiently close Xi n and Xf. For more general considerations 
we divide the trajectory into small parts of length Ax for which the above 
expression is still true and then perform a summation (or rather an integra- 
tion in the limit Ax — > 0) over all possible trajectories. 

In more explicit terms, what we do is to express the propagator kernel 
(and the Green function) by such path summation. As we know from Huy- 
gens principle, the propagator responsible for the evolution of ^(qi n ,Un) ~^ 
i/j(qf,tf) can be expressed as a composition of infinitesimal propagators 

. JV-l 
K{q in -> q f , t in -*tf)= / dqidq 2 . . . dq N _i JJ K(q n -> q n+1 , t n -> t n+1 ) 

J n=0 

with g = q%ni qN = qf- Since, as we already know, for small intervals 

K(q n -> q n+ i, t n -> t n+1 ) = (q n +i\exp[- ' lH ^ t ]\q n } 

separating the kinetic and potential term in H and neglecting 0(At 2 ) terms, 
we may approximate each of the infinitesimal propagators as At — > with a 
factor in the form 



/ m 



2nihAt 



-1/2 



exp 



iAt ( m ( Ag x 



h \ 2 \ At 



- V 



22 at this stage we are already dealing with a semi-classical method since points in purely 
quantum phase space are ill-defined due to Heisenberg uncertainty principle 



30 



which gives for the overall kernel 



K(q 



in 




J "/ \2mhAt) 



(27) 



n-i [i r l f 

x Yl dq j (2nhAty 1/2 exp - L(q,q,t)dt 




Note that most 'weighted' paths are namely the classical trajectories, for 
which the action functional is extremal. On the other hand, a path, partici- 
pating in the above distribution does not necessarily have to be smooth, as 
can be seen in the example of Brownian motion. 

From path integrals one may readily retrieve information about transition 
coefficients, that are propagator entries, but also, via Fourier transform one 
may recover the Green function. Calculation of path integrals in practice , 
however, is not that simple. In the trivial case of quadratic Lagrangian, the 
expression is reduced to the standard Gaussian integral, that is computed ex- 
actly. For the more generic case, however, we need to rely on approximation 
procedures, such as the stationary phase approximation, or the perturbation 
series expansion, both considered in the following. 

3.5 Perturbation Theory 

There are but few known problems in quantum mechanics that allow for ex- 
act solutions and for practical matters, we need to be equipped with good 
approximation methods. Perhaps the most used one (not only in quantum 
mechanics, but in classical dynamics and QFT as well) is the perturbation 
method. 

The simplest illustration of the method^) starts with an unperturbed Hamil- 
tonian H , for which the stationary Schrodinger equation is solvable by an- 
alytical methods, or for some reason we know the spectral decomposition 



Then imagine that we want to go a bit farther and study the eigenvalue 
problem of the operator H = Ho + V, where the term V is small compared 

23 this is the easiest case with time-independent hamiltonian with discrete spectrum, but 
also the most relevant to the topics discussed here 
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to Hq in the sense that \Vki\ -C \E% — Ef\ where the matrix coefficients are 
computed in the initial basis {ip k }- This restriction will be justified in the 
following. 

Now, as usual, one expands the unknown eigenfunctions ipk in the initial 
basis as 

3 

and substitutes into the equation for H. What one gets is 

k m 

Multiplying both sides by ip® and integrating yields 

(JB-£?)c, = (28) 
i 

which is a homogenous system for the Cj's which leads to an infinite series of 
relations for the perturbed energy states Ej in terms of the unperturbed ones 
E% and the matrix entries of V. However, considering the fact that the latter 
were chosen to be small, we may expand E k in a series E k = E®+E' k +E k +. . . , 
where the j th term is a polynomial of order j with respect to the V^'s and 



cut that series at an appropriate place 24 !. The same procedure we do with 
ip, or rather, with the cVs. Then we restrict to the eigen-state ijj k by taking 
c° k = 1 and c° = for j ^ k in the above equation. Substituting all this 
in the equation we may immediately determine the first correction E' k . For 
k = i it gives 

E\ = V kk (29) 



and for i ^ k, obtain the first correction to the eigenfunction 



(30) 



In order to have this correction small, we really need the estimate we imposed 
in the beginning. Then the perturbed egen-state ip is given, up to a first-order 
corrections, by 

^~^° + £l^^ (31) 
i^k ^ 



Hn practice, it usually suffices to take the first two or three terms 
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Note that the correction term lies in the orthogonal complement of if)^, so 
the new state ipk is normalized up to second order corrections. 
Following the same procedure, we easily obtain the next few terms 



E 't = £ 



IV, 



ik 



1 2 



' & h < B ° - E -°»< B ° - E °'> ~ u ~ E - )2 

The results obtained here can be generalized to the case when H has a con- 
tinuous part of the spectrum as well. Then summation is simply replaced by 
integration. 

In the case when the spectrum of Hq is degenerate we still have to modify 
our procedure a little, due to the arbitrariness of the states in the degenerate 
subspace. It happens, that this arbitrariness is eliminated once we impose 
the condition that states should be altered only by small variation. Such 
'correct' states are defined as combinations of the initial ones with coeffi- 
cients c° , . . . c° where s denotes the degree of degeneracy of the eigenvalue 
Ef. With this substitution, equation (128]) . restricted on the corresponding 
invariant subspace gives 

J2(Vu k -E'5 Uk )cl=0, t,k = l...s (33) 

k 

which as a linear homogenous system for the c?_'s has nontrivial solution if 
and only if 

det(V iik -E'8 iik ) =0 (34) 

which on its side is a polynomial of degree s in E' and its s roots give 
the energy shifts within the degenerate state up to first order. In order to 
determine the 'correct' eigenfunctions, we substitute these roots back in the 
matrix equation in order to get all the c° fc 's. Then, for the perturbation of 
the eigenfunctions we find 

V ikl V u 



lk V u -V ikik ^E°-El 



Evi k ivii . , . 



It may happen that some of the corrections (or even all) are zero. Then one 
proceeds with higher order perturbations. On the other hand, if the roots 
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are distinct, so are the perturbed levels and this is how perturbation often 
cancels degeneracy. 



Now, for the sake of consistency, we need to consider time-dependent case as 
well. Let 

in the first place and H = H + V(t). Now we expand, just as before, the 
corresponding solution in the form 

k 

and substituting into the equation for H, we end up with a ODE for the 
coefficients a,ki(t), namely 

ihau = ^2 v km(t)a ml (35) 

m 

Then, proceeding as before, we get for the first correction to the I th state 

iha'u = V k i{t) 

or 



»«(*) = ^ / V M (t) dt 



There are numerous interesting applications of perturbation theory in diverse 
branches of classical and quantum Physics. We only point out here one more 
important aspect that goes far beyond the basics of quantum mechanics, 
considered here. 

Let us take the kernel of the propagator for the transition q in — > qf of a free 
particle and denote it by K (q in — > qf). We already know that this kernel 
can be represented as path integral in the form 

K (q in -> q f ) = J exp ^ J £ (q, q) dtjVq(t) 

where Co(q, q) = mq/2 is the classical free-particle Lagrangian and V denotes 

the measure of the path integral, as found in the previous section. 

Now, let us perturb this Lagrangian with a small potential term £(q, q, t) = 
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£o — V(q, t), V <C C and pose the question of finding a reasonable approxi- 
mation of the perturbed propagator kernel Ky. 

Assuming that V is small compared to h in the whole time interval, we may 
expand the ^/-dependent exponent factor in the path integral in Taylor series 

exp (l I ' V(q,t) dt \ = 1 + % - ! ' V{q,t) dt - ^ ( V V{q,t) dt\ +... 

\ J ti n ' J tin \" tin ' 

which leads to an expansion for the overall kernel in the form 

K = K° + K' + K" + ... 

For the calculation of the distinct terms it is helpful to note that Ky is a 
fundamental solution of the Schrodinger equation in the sense that 

/ d h 2 \ 

y h df f + 2^ V "f ~ Ky = in5 ( t f- t in)^f - Qin) (36) 

and so is K° with V = 0. This allows for applying the parametrix method, 
briefly discussed below in the case of the heat equation, to expand Ky as an 
infinite convolution product geometric series^], which allows for writing the 
n th correction term in the form 



i rtf rt n -! i>gi 

K {n \q m ^q f ) = — / . . . / if°(g n -l - 9/M<2n-l, t n -i] 

\ lrL ) Jt n -i Jt n - 2 Jq in 



X 



-1 J t n -2 

xK°(q n _ 2 -> q n -i) V(qi,ti)K°(q in q t ) dt n dt n -\ ...dt x 

This term is interpreted as a propagator of a particle that has been scattered 
n times by the potential V within the interval (ti n ,tf). This interpretation 
is illustrated by the so called Feynman graphs, that are used as a means of 
calculation as well. 

The zeroth term in particular describes the probability that the particle 
remains free without noticing the presence of potential. Then K' v stands for 
a one-collision process etc. Resuming, if possible , all (or at least a sufficient 
number) of the above terms, we get the demanded propagator Ky (or a good 
approximation of it). Such summation leads to the integral equation for Ky 

1 f tf 

Ky(q in -> q f ) = K°{q in q f ) + — / K°(q -> q f )V(q, r)K v (q in -> q) dr 

J tin 

(37) 

which can be obtained also directly, using the differential equations for the 
two kernels. 



2r, 



known as Voltera 
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3.6 Scattering Problems 



We end this section with a short discussion on the stationary Schrodinger 
equation defined on the line R, which for simplicity we write in the form 

(H - k 2 )ip = (38) 

d 2 

where H = — — - + V(x) represents the re-scaled hamiltonian with a finite 

ax 2 

potential (or 'scattering') term V and the spectral parameter k 2 stands for 
the energy. 

The whole space of solution can be expanded over a basis, chosen in the 
following way: consider a solution %p r + with that tends to e lkx as x — > oo 
and another such i]f_ which tends to e~ tkx . These two independent solutions, 
sometimes called right Jost solutions, represent, in physical terms, two plain 
waves traveling to and from plus infinity, respectively. There is, however, 
another option to construct such a basis by switching from left to right infinity 
ip l ± ~ e ±ikx as x — > ±oo. Since the full space of solutions of the above second- 
order ODE is two-dimensional, these two sets must be linearly dependent. 
We shall construct such dependence in an instant. Consider for example a 
solution that behaves like an outgoing plain wave as x tends to infinity and in 
the other limit, it should be given by a linear combination of the two 'basic' 
states, or 

i> ~ a(k)e lkx + b(k)e~ ikx 

where a(k) and b(k) are holomorphic in C/{0}. The physical picture repre- 
sents two asymptotically plane waves - one initial and one that is 'scattered 
backward' by V(x). Alternatively, we may take the initial wave traveling 
from the right to the left - a solution that tends to e~ lkx as x — > — oo and 
as x — > oo fulfills the limit 

~ p(k)e lkx + a(k)e- lkx 

The same considerations hold here as well. Moreover, both <fi and ip satisfy 
the property that the Wrontzkian of the function and its complex conjugate 
is constant. For the proof we only use that both functions are solutions of 
the stationary Schrodinger equation, given above. Take for instance ip and 
substitute then multiply by ip. The same procedure we repeat with inverting 
ip and ip - note that ip is still a solution, since the Schrodinger operator is 
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self- adjoint. Then we subtract the two terms and finally obtain 




which means that 



W(ip,ijj) := ( 4>-rip — i^-rii ) = const 
V ax ax J 




The same certainly holds for as well. 
From this single relation we get 




On the other hand, in the same way as before, it is obvious that W((j),ip) = 
const, which immediately leads to a = a and hence |6| 2 = |/3| 2 = l + |a| 2 . 
Now let us think in more physical terms and regard the whole problem from 
the initial wave's point of view - this means, above all, choosing it's amplitude 
to be unity or, in other words, simply dividing by the non-zero factor a(k), 

b 3 
so that the reflected term adopts a factor r + = — ( respectively r_ = — ), and 



the transited comes with t = —. The letter are called reflection and transition 

a 

coefficients respectively. Now we define the scattering matrix of the problem 
as 



It clearly gives the scattered state (in both directions) expanded over the 
basis of Jost solutions. Moreover, the matrix S(k) is unitary - that is un- 
derstood most easily if we look at it as representing an evolution operator, 
acting on quantum-mechanical states. It also appears to be real in the sense 
that S(—k) = S(k). 

There is one more significant operator involved in the very root of our prob- 
lem. Essentially, this is given by the monodromy matrix A4, defined in 
the following way. We start with a solution ip that has the form of a single 
plain wave as x tends to minus infinity. Let ip Q be the unique solution of 
the free-particle (Helmholtz) equation, that coincides with if) at this limit. 

26 one may check that their zeroes are purely imaginary 



a 



a 



1 




37 



Alternatively, we denote by ipQ the solution that coincides with ip in the limit 
x — > oo. Then the matrix M. is defined as the unique linear isomorphism 
which satisfies 



In order to be more explicit, we express a generic monodromy matrix by 
means of the transition and reflection coefficients. We have 



This turns out to be a SU(1, 1) matrix, as can be derived from a simple 
consideration: construct a real basis (involving sine and cosine instead of the 
complex exponents) in the space of solutions of the free-particle equation. 
Then, the monodromy matrix (in this basis) should be area-preserving^!, 
and hence, SL(2,M), which in the complex basis, as we know, takes the form 
of a SU(1, 1) - matrix. This is enough for us to prove the above formula. We 
just note that it transforms the plain wave e tkx + r + e~ tkx into the plain wave 
te lkx which immediately gives the first column of the contra-gradient matrix, 
transforming the basis e± = e . For the second one we only have to take 
the complex conjugate to obtain finally 



then we easily obtain the monodromy matrix, exploiting the fact that to 
invert a SU(1,1) matrix, we only have to exchange the diagonal elements 
and take the non-diagonal with opposite sign. 

With this formula in mind it is easy to see that the sum of the transition 
and the reflection amplitudes is always unitj@, which is the fundamental 
principle of matter preserving in physics. 

To sum up, we built in this paragraph a tool through which to each finite 
scattering potential one may assign a monodromy matrix in the group of 
isometries of the unit disk (or, equivalently, the upper-half plane, depending 
on the basis chosen). Alternatively, one may assign a SU(2) - scattering 
matrix, which in many situations possesses additional symmetries, giving 

27 that is because it can be thought of as being adjoint with the phase-area preserving 
flow, acting in the space of 1-jets of solutions of the initial equation 

28 the amplitudes are defined as the squared moduli of the coefficients and have the 
physical meaning of probabilistic densities for transition and reflection, respectively 



Vo = M1>\ 
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birth to a means of reducing interaction problems - a common practice in 
Quantum field theory for example. 

Sometimes it is sensible to ask weather one can restore the potential from 
the scattering data: the S'-matrix, together with the finite spectrum (the 
finite set of negative eigen-energies, that correspond to bound states) and the 
divergent points of the transition coefficient. This is the inverse scattering 
problem that turns out to be solvable in some cases and in addition gives 
rise to a powerful tool of resolving a special kind of non-linear differential 
equations, called for historical reasons solitons. 
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4 Towards the Theory of Integrable Systems 



The theory of non-linear integrable equations (or infinite-dimensional hamil- 
tonian systems) is relatively new, but has already evolved far and spreads 
over various branches of modern mathematics, such as group theory, func- 
tional analysis and others. Our aim here is only to outline the close relation 
between this topic and our main concern - spectral theory. 



4.1 Hamiltonian Formulation 

Our base point so far was a quadratic hamiltonian (or lagrangian) description 
of a classical dynamical system, naturally leading to linear ordinary differ- 
ential equations. The problem of integrability from standard ODE theory 
was reduced to finding a full set of functionally independent and pairwise 
commuting integrals of motion. Conceptually the same thing refers to the 
non-linear case, but this time the action-angle basis is infinite-dimensional 
and the problem of finding one becomes much more complicated. 
Our only constructive example until the end of this section will be one of 
the first well examined non-linear PDE's - the Korteweg - de Vries (KdV) 
equation 

u = 6uu' - u'" (39) 

where ' stands for time derivative, while ' means differentiation with respect 
to x. 

First of all, let us introduce our algebra of observables as the space of real- 
analytic functionals F[u] and then introduce on it a Poisson structure, defined 
by the bracket 

W = (40) 

where the <5-sign stands for functional derivatives. 

It is clear that this bracket defines also a symplectic structure in the usual way 
and we can easily build a Hamiltonian formalism with its help. In particular, 
once we have chosen a Hamiltonian functional F[u], we recover the evolution 
equation of an observable via the Poisson bracket in the standard form 

f = {f,m (4i) 
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If we choose the hamiltonian in the form 



dx 



(42) 



for example, we recover the KdV equation 







(3u 2 - u") 




where the last term is obtained via integration by parts. 
With the help of the Poisson structure, one generalizes the Liuville-Arnold 
theorem to the infinite-dimensional case and claims that a system, defined 
by the hamiltonian H[u] is integrable if and only if it allows an infinite set of 
functionally independent and pair-wise commuting integrals of motion I k [u] , 
that are functionals, commuting with the hamiltonian in the Lie algebra, 
defined by the above bracket. 

In our example we have first Iq = H, the second integral Ii[u] = J u(x) dx is 
concerned with the structure of the Poisson bracket and more specifically, the 
presence of a spacial derivative in it. The third preserved quantity follows 
from translation invariance of the hamiltonian functional. On^ immediate 

consequence of this invariance is the vanishing of the term / ———dx that 



is the first correction for the small translation parameter expansion of H[u]. 
Now it is obvious that if we define 



As we shall see shortly these are really just the first few representatives of an 
infinite set of integrals of motion, each of which can be itself thought of as a 
hamiltonian of an integrable system. Therefore these system come in infinite 
hierarchies defined by the set of hamiltonians Ik[u]. For example, from the 
above we easily see that in the hierarchy of the KdV equation we have also 
the evolution equations u — u', u — 2uu' etc. 

4.2 Lax Representation and Iso-spectral Deformations 

Using again the KdV equation, we construct an equivalent operator repre- 
sentation, known as the Lax representation. To begin with, we consider the 



R 



Su dx 
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Schodinger operator for the potential u from the KdV equation 

L = -d 2 x + u(x,t) 

that is a special case of Lax operator we are going to investigate here and 
consider also the equation 

L = [L, M) (43) 
where the operator M is given by 

M = A&1- 3(ud x + d x u) (44) 

A simple calculation shows that the left hand side multiplies a test function 
by it, while the remaining factor from the right-hand side is exactly 6uu' — u'" 
which recovers the initial partial differential equation for u. 
Now let us see what else we can retrieve from it. First of all, consider a 
solution <f) of the scattering equation 

L(f) = k 2 (f) 

with the asymptotic of a left-going wave at minus infinity Next we differ- 
entiate the latter with respect to t. Substituting into the Lax equation we 
obtair@ 

(L- k 2 )(0- M<j>) = 

so the function = — M0 is also a solution, corresponding to the same 
eigenvalue k 2 . The asymptotic behavior of completely determines that of 
0. We only substitute to see that 

~ Aik 3 e~ lkx as x —> oo 

Since this behavior does not depend on t, substituting the explicit expression 
for gives the time- dependence (|47p of exactly as promised in the begin- 
ning. 

Let us now go back to the scattering problem for the Schrodinger equation 
from the previous section , where the quantum hamiltonian coincides with 
the Lax operator. Note that we let this time the potential depend on an 
additional 'time' parameter t, so the solutions and all scattering data will 

29 here we use the fact that the spectrum of L is preserved with time that is to be 
explained below 
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depend on t as well. It turns out that in some cases, a quite complicated 
non-linear evolution of u(x,t) may result in a linear transformation of the 
^-matrix entries. Consider again the KdV equation 

u = 6uu' — u" (45) 

This quite cumbersome evolution equation for the potential naturally induces 
a quite simple transformation rule for the scattering coefficients, namely 

a(k, t) = 0, b(k, t) = 8ik 3 b(k, t) (46) 

and the 'scattered wave' evolves like 

ijj = (Aik 3 - M)V> (47) 

With this in mind, it is quite natural to use the correspondence potential- 
scattering data and transform our non-liner problem into a linear one. 
In order to see how the 5-matrix evolves with time, all that one needs to do 
is to substitute in (T4"T|) the right-infinity limit for (f>: 

(j)(x, k, t) ~ a(k, t)e lkx + b{k, t)e~ ikx 

All these equations are linear and easy to solve. In order to benefit from this, 
however, one needs to translate backwards the information - from evolved 
scattering data, to potential at time t. This leads to the Gelfand-Levitan- 
Marchenko integral equation of the form 

POO 

K(x,y) + F(x + y)+ K(x, z)F(z + y)dz = (48) 

Jo 

where the function F is defined by the formula 

F (x) = Y — ( r + — / r(k)e lkx dk 

All the functions on the right hand side are completely determined by the 
scattering data but we are not going into detail with this dependence. 
Instead we refer to [29] and especially for the inverse scattering method 



:S0 



including the spectrum, the iS-matrix and the scattering poles residua 
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[To] . Then the solution of the above integral equation recovers the potential 
function by the simple formula 

u(x) = -2K'(x,x) (49) 

One immediate consequence of the the skew-symmetric nature of M is the 
unitary evolution for the Lax operator, which means that at time t, we have 

L{t) = U{t)L{Q)U\t) 

where the unitary operator U has the form U(t) = exp(-Mt) and inversely, 
M = UU^ 1 . This means that L(0) and L(t) are unitary equivalent and there- 
fore the spectrum together with all spectral invariants is preserved with time. 
This observation has two straightforward consequences. One of them is the 
use of the spectral invariants as generators of infinitely many pair-wise com- 
muting integrals of motion that we need to resolve the system in Hamiltonian 
terms. For the special case of the KdV equation this procedure has a lot to 
do with our further considerations on spectral theory. For the time being 
we only mention that first few integrals of motion for the KdV equation are 
polynomials of u and its spacial derivatives, given by 



1 1 = — J udx 

u 2 dx 



u' 2 + 2m 3 ) dx 



ir2 



5u 2 u" + 5m 4 ) dx 



The second consequence we would like to mention here is the following. Con- 
sider the geometry of a smooth manifold with a Laplace operator L defined on 
it (u might represent curvature term for example, and the ordinary deriva- 
tives might be replaced with covariant ones). Then we study all possible 
smooth deformations on the manifold and the way they reflect the operator 
L. We call one such deformation iso-spectral if it preserves the spectrum of 
L. In this context, each such deformation generates an integrable partial dif- 
ferential equation (infinite-dimensional hamiltonian system) and vise versa. 
This dualism between iso-spectral deformations and integrable system may 
be used in several directions but this investigation, as they often say, goes 
far beyond the scope of the present work. 
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5 Calculus 



After providing some physical and geometrical insight, it is time to introduce 
some of the the basic tools and results used for practical calculation. 

5.1 Laplacian Spectra 

The use of spectral theory in physics is so common that we often put a sign of 
equality between certain mathematical spectral problem and a whole branch 
of the physical theory as it is in the case with the Maxwell, Schrodinger, 
Klein-Gordon, Dirac, Poisson equations etc. 

In this section we give the mathematical roots of the beautiful physical the- 
ories derived by the equations, mentioned above. In order to do this, we 
introduce the most typical classes of PDE's, referred to as the equations of 
mathematical physics, and give a first hint how their spectral properties are 
linked to the geometry of the underlying manifold. 

This is going to become clear from the examples that follow. We along with 
them give some facts from the theory of partial differential operators. 
To begin with, each such operator (of order k), acting on a manifold X has 
the formal 

P(x,D) = J^)D a (50) 

\a\<k 

and can be represented in terms of its symbol 

a{P){i) x = P{x,0 := J «W* (51) 

\a\<k 

This is simply a polynomial in the dual variable £ = ^dx 1 G T*X, obtained 
by a Fourier transform. Its homogenous part of highest degree 

a L {P) := 

\a\=k 

is called the leading symbol of P and defines an invariant^ homogeneous map 
of degree k from T*X to the space of coefficients (in the case of matrix- values 

31 here and below a is a multi-index, meaning that D a — d^d^ 2 . . . d" n (m is called 

n 

'the length' of the multi-index and |a| — ^""^ otj - its 'volume'). 
32 in fact it is the only covariant part of the symbol 
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coefficients, this is the space of homomorphisms between the corresponding 
vector bundles over X). 

Now P is said to be of Laplace type if there exists a metric g on X such, that 
the leading part of the symbol can be represented as a scalar square with 
respect to this metric ctl(P) = \ \C\\g- 
P is said to be of Dirac type if P 2 is of Laplace type. 

There are numerous examples of such operators in almost every branch of 
modern mathematics. We mention here some of the simplest and most intu- 
itive ones. 

1. Consider a smooth compact and orientable manifold M without bound- 
ary and regard the spaces of differential k- forms Q k (M) = A k (T*M) 
and the differential, acting on them 

d : Q k {M) -> Q k+1 (M) (52) 

The property d 2 = turns (|52|) into a differential complex with co- 
homological groups^ 

H k (M,R) := ker(4)/im(4-i) (53) 

the so-called de Rham co-homologies. In the spaces Q k we have also 
a well-defined scalar product, thanks to the Hodge duality * : Q k — ► 
Qm-k^ m _ ^m(M). Namely for two fc-forms u and <fi we have the 
coupling 

{uj,4>) := / u A *0 (54) 
Jm 

which allows for defining the conjugate differential (or co- differential) 
5 : Q k —> Q^ 1 by the equality (w,6<f>) = (du,</>). Then, since 5 is 
immediately nilpotent 5 2 = 0, the operator 

A H := d5 + 5d 

is by construction self-adjoint, and its leading symbol is given by a 
scalar square with respect to the above metric (this is easy to verify 
using a suitable basis). 

Ah '■ Q k — » Q k is called the Hodge Laplacian and the restriction of its 
kernel to the space of fc-forms has the dimension of H k (M, R) - there 

33 <4 is in this case the restriction of d to Sl k (M) 
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is exactly one representative in each cohomological class. This is the 
statement of the famous Hodge - de Rham theorem, providing the first 
example of the great topological and geometrical significance of the 
operators we are concerned with herd^l. 

2. Consider now a (pseudo-)Riemannian manifold with Levi-Civita con- 
nection V. It defines the Laplace- Beltrami operator, by the formula 

A LB = V l V, = -j^d, [VW\9 ij dj) (55) 



\9 

where are the entries of the metric tensor and \g\ - the modulus of its 
determinant. This operator happens to coincide with the one from the 
previous example only when acting on scalar functions, or in the trivial 
case of flat manifolds. For 1-forms, they disagree with a term, equal to 
the scalar curvature of the manifold and for higher rank tensors, this 
relation becomes more complicated. 

For physical applications it is often convenient to use the fact that a 
generic Laplace operator 

P ■= -g(d,d) oid + a k d k + b 

can be represented in the form 

P:=-g{D,D)-E (56) 

where D and E are a Riemannian connection and an endomorphism, 
chosen as follows: 

A* = \g* («* + 0"ry 

E = b- g ij {did* + AiAj - A k T) k ) (57) 

The above note makes it now quite obvious why in classical field theo- 
ries we talk about curvatures and and potentials as if they are equiva- 
lent. 



in electrodynamics and fluid mechanics the operators d and 5 play central role and are 
usually denoted by rot and div respectively - more precisely, they act in this way on one 
forms, while on m — 1-forms it is the other way round, since 5 = ± * d* and the Hodge 
duality on forms holds 
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Exploiting some tricks from differential geometry we obtain the promised 
relation between the two Laplace operators defined so far, at least for 
1-forms and vector fields. Namely, in this case we have 

= (d6 + 5d) ]nl = -V*Vi + Ric (58) 

In fact this is only the first of a series of formulas, relating the Hodge 
laplacian Ah with the Laplace- Beltrami operator Alb, known as Bochner 
- Weitzenboch formulas. In a flat space (possibly compactified) cer- 
tainly both laplacians coincide. 
One may go a bit farther and prove that 

iA|£| 2 =|V£| 2 ±ffic(£,0 

in which the 'plus' sign refers to harmonic vector fields £ for which the 
tensor V£ is symmetric (they are locally given by gradients of smooth 
functions), whereas the 'minus' sign is taken on the space of Killing 
vector fields - the ones for which V£ is skew-symmetrid 35 !. 
Integration over M in both cases gives zero and when in the right hand 
side we have two non-negative summands, they must both vanish inde- 
pendently. In this way we easily conclude that if M is closed, and the 
Ricci tensor is non-negative, then all harmonic one-forms (and respec- 
tively vector fields) are paralleS In addition to that, if Ric is positive 
definite at least at one point (and still non-negative everywhere else), 
then £ is bound to vanish there due to non-degeneracy, but since it 
is parallel, this means that it vanishes everywhere. In this case there 
are no non-zero harmonic one-forms or vector fields on M at all. This 
result is known as the Bochner vanishing theorem. 
These quite simple considerations provide some useful information con- 
cerning the topology of the underlying manifold. First of all, if \ \Ric\ \ > 
then the dimension of the maximal flat subspace imbedded in TM (in 
the case of symmetric spaces, this is their rank) equals the first Betti 
numbeiffi Pi = dim i7 1 (M,M). Therefore we have (3\ = fi m -\ < m — 

35 recall that a Killing field X is by definition a field, generating isometry, and satisfies 
V V X^ + V^X U = 

36 we call a £ f2 1 (M) parallel with respect to the connection V if Vx(cr) =0VXe TM 
37 since the above considerations about vector fields naturally refer to one-forms as well 
and one may easily verify that all parallel sections are harmonic 
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dimM. Equality holds only for tori, while on spheres we obviously 
have no global parallel sections. 

In the case of strictly negative scalar curvature, the above equation 
could be used to verify that there are no Killing vector fields on M. 
Therefore, the group of isometries could only discrete and thus, due to 
compactness, finite. This result may be used directly to prove that all 
geodesies in this case are isolated. Another possible use wold be in the 
attempt to obtain hyperbolic manifolds as ball quotients. 

3. Now suppose that we have the flat Minkowski space IR^ 1,3 ) and the 
D'Alambert wave operator 

= 81-82-83-84 

acting on it. This is a well defined Laplace type operator that possesses 
a square root. Namely, we may construct a Dirac operator with the 
help of the matrices^! 7°, satisfying the relations 

YY + YY = -2^id (59) 

Here g is the flat hyperbolic (Lorenzian) metric and our operator is 
simply given by 

p = 7 ^ + 7° 

Note that the existence of Laplace operators is guaranteed by the metric 
and the smoothness of the manifold, whereas their square root might 
not exist. Speaking in general, the presence of a Dirac operator, acting 
on X is a matter of whether the vector bundle of interest admits a 
spin stricture, or a representation of the Clifford algebra, which is the 
universal unital algebra, generated by vectors in M n subject to the 
relations fl59|) . If so, we have a double cover of groups: 

Z 2 -> Spin(n) -> SO(n, R) (60) 

The Spinor groups in the low dimensions happen to coincide with the 
universal covering groups of the special orthogonal groups. For example 
Spin(3) = SU(2) = § 3 , or Spin(4) = § 3 x§ 3 etc. In order to investigate 
the problem of existence of spin structures on manifolds, one needs to 

8 usually referred to as Dirac matrices 
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study some topological properties of the corresponding principle bundle 
with the help of characteristic classes. Namely, we use the notion of 
Stief el- Whitney classes oJiiV) G H l (X, Z 2 ), where V is a real vector 
bundle over X. For example, if io\ 7^ (and certainly H 1 (X, Z 2 ) 7^ 
0), then the bundle V is not orientable, as we know from elementary 
algebraic topology In the case it is, the second class u 2 is an obstruction 
for the existence of spin structure. Furthermore, the Abelian group 
H 1 (X, Z 2 ) parametrizes the inequivalent spin structures on V. 
Using this apparatus, we easily obtain that the spheres of all dimensions 
are spin as well as odd-dimensional complex projective spaces. In the 
real case MP m is spin if and only if m = 4k — 1 and has two inequivalent 
spin structures as H 1 (RF 4k -\Z 2 ) = Z 2 . For a spin manifold we have 
at least one Dirac operator and a spinor connection, that agrees with 
it. In our case it is given by D = 7^V y . Then D 2 is known as the spin 
Laplacian and might be showrj^l to be of the form 

A spm := Tr(V 2 ) + ~ (61) 

with k(x) being the scalar curvature of the underlying manifold. This 
formula helps for retrieving results analogous to these from the previ- 
ous example. Namely, we have that if the scalar curvature on M is 
non-negative, all harmonic spinors are parallel and if at least at one 
point it is positive, then there are no harmonic spinors at all. 

The three major classes of second-order linear PDE's, concerned in mathe- 
matical physics in which Laplace type operators are involved are: 
elliptic equations of the type A<p = p(x), where the Laplacial is taken with 
respect to a metric with positive (elliptic) signature. The function p is usually 
interpreted as a volume distribution of the charge and <f>(x) is the potential 
field. Sometimes we regard the homogeneous equation with p = to investi- 
gate stationary processes (electrostatics, distribution of heat etc.) 
parabolic equations involve more (for example time) derivatives of lower or- 
der, which are absent in the leading symbol. Therefore, the heat equation for 
instance, known also as continuity or diffusion (A + d t ) <p — is considered 
to be with degenerate leading symbol. 

The last class is the class of hyperbolic equations of the above type □(/> = p(x) , 
39 again using Bochncr technique - the above formula is due to Lichnerovicz 
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but this time with hyperbolic metric. One of the most popular examples is 
the D'Alambert wave equation, considered in the previous section. In par- 
ticular problems the save equations is reduced to a spectral equation of the 
type A0 = k 2 (p, known as the Helmholtz equation. 

5.1.1 Boundary Value Problems 

In the case of closed manifolds the Laplacian appeared to be self- adjoint by 
construction, but this property is no longer guaranteed when introducing a 
boundary. According to the Green formula, one has to encounter boundary 
terms as well. In order to get rid of them and to keep the symmetry of A, 
one imposes appropriate conditions for the behavior of the solution and/or 
its normal derivatives at the boundary. 

As we may see from the following considerations, there are some naturally 
preferable ways to do this. 

The Green Function and Layer Potentials For the so called equations 
of mathematical physics, it is possible to reduce the problem to the bound- 
ary, using single and double layer potentials, defined by means of the Green 
function and thus transform the initial PDE into an integral equation on 
the boundary. 

Let L(x) be a general self-adjoint PDO with varying coefficients and A(x,£) 
represent its integral kernel, given by the Fourier inverse of the characteristic 
polynomial. Thus we have 



The kernel G(x,£), referred to as the Green function of L(x) is determined 
by the sole definition of inverse L o L^ 1 = id, which leads to 




for any solution ip and we may write formally 



rj) = L~ l F(x) 



where the formal inverse is an integral operator. We write 




L(x)G(x,0=6(x-0 
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Let us now take L to be the Laplace-Beltrami operator V 2 in a bounded re- 
gion Q, and apply the Green theorem for the Green function and an arbitrary 
solution ip. We have 



Jan 



du du 



da 



where d/du = zaV is the inward-pointing normal derivative at the boundary 
Considering that ip is a solution and G is the kernel of the formal inverse, we 
obtain from the above 



G.Fduj + 



an 



du du 



da 



(62) 



In order to solve this integral equation, we need to get rid of some of the 
terms. This is being done by imposing appropriate boundary conditions. 
There are three types of boundary conditions used in this situation: 

1. Dirichlet boundary conditions demand the restriction of ip on dfl to 
be equal to some known function /. The solution is obtained if we set 
G = on the boundary in order to eliminate the term involving the 
normal derivative of tp on <9f2. We finally have 



^(0= / G.Fduo- I 
Jn Jo: 



an du 



(63) 



2. Neumann boundary data prescribes the behavior of the normal deriva- 
tive of the solution at the boundary d u ip\QQ = f. By analogy with the 
previous case, one may attempt to set d u G = on dfl, but such so- 
lution would be ill-posed, as one may see, substituting into the Green 
formula. We may easily get rid of this inconvenience, and still have 
d u G\Qn = 0, if we let the Green operator to be inverse just modulo a 
constant term 

V 2 G = 5 + C 

The constant we really need is the inverse volume of Q (taken with 
'minus' sign). The Neumann BVP then reduces to 



=< ^ > + I G.Fdu - I G.fda 
Jn Jan 

with < ip >= Vol(fi)" 1 J Q ip du 



(64) 
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3. The Robin cas cj is obtained if we set a combination of the solution and 
its normal derivative at the boundary to be equal to a known function 
d u ip + Sip = f on dQ. It is reasonable to ask the Green function to 
satisfy d u G + SG = at the boundary. Then, applying directly the 
Green formula, one gets again 

V>(£)= / G.Fdu- I G.fda (65) 
Jq Jan 

The boundary terms, appearing in the above expressions are usually referred 
to as layer potentials. More precisely, the single-, respectively double-layer 
potentials are defined in terms of the integral operators 

S£(k)[f](x) = I G(k,x,y)f(y)dy((P) 
Jan 

r bg 

Dl(k)[f](x) = & —(k,x,y)f(y)dy( ( p) (66) 
Jan ° v 

Now the question of how we build Green functions, satisfying given boundary 
conditions, naturally comes to mind. To understand this, one first takes a 
more general fundamental solution, known as the free Green function Gq and 
then adds correcting function in order to fulfill the demanded behavior at the 
boundary. 

Note that the non-homogeneous boundary conditions do not constitute a 
vector space - for example a sum of two such solutions is no longer a solution, 
so one may think that he has come upon an obstacle. However it is easy to see 
that all non-homogenous solutions of a given equation of the above type may 
be obtained form the homogenous ones via modification of Green function, as 
we already saw, so that it suffices to consider only the homogeneous casa 41 !. 
Let us come back to our initial example - the Laplace-Beltrami operator, 
acting in a bounded domain. In the two-dimensional case, the free Green 
function is known to be proportional to lnr, r = a/x 2 + y 2 , otherwise Go ~ 
r 2 ~ d , d ^ 2 - for each dimension there is a characteristic divergency, which is 
rather topological in nature, but easy to find by analytical methods as well. 

40 called also 'mixed', or 'modified Neumann' 

41 the same refers to the question of homogeneity of the equation itself - if ^ is a solution 
of the homogeneous equation, then ip = (f> * / satisfies non-homogeneous equation with / 
instead of zero on the right, since the convolution product is compatible with the operators 
we are concerned with here 



53 



It is obvious that these Go's are harmonic, except for the singularity at the 
origin. Adding a regular harmonic function h does not affect that property 
and that is exactly how boundary conditions are obtained. 
For the heat equations, it is easy to show via Fourier transform that the 
fundamental solution has kernel in the form 

G (x,y,t) = (Ant)- d / 2 e- lj ^ jil 

A particular interest for us in the following will be the two-dimensional 
Helmholtz equation V 2 ip + k 2 ip = 0. Its free Green function is shown to 
be given by a Hankel function of the first kind 

G (x,x') = ^H^lkix-x 1 )} 

5.1.2 Estimates and Properties of the Spectrum 

Here we study some of the properties that a Laplace spectrum on a compact 
Riemannian manifold should possess. Most of them could be found either in 
[Ml or in [21]. 

First of all, let us consider again the eigenvalue problem 

A0 = X k (f) k (67) 

with the corresponding boundary condition (D,N or R). Then we have for 
the spectrum growth rate the following estimate: 

1. > 3n(e) E N, n™~ e < \ n < n™ +£ , where m is the dimension 
of the underlying manifold M 

2. let Cfc denote the Fourier coefficients in the expansion of any L 2 - func- 
tion ip over the basis (f>k of eigenfunctions of A. Then ip is smooth if 
and only if lim n k c n = 0, VA; G Z 

3. if in M we have a Riemannian connection and \\ip\\k denotes the sup- 
norm of the k~ th covariant derivative of the L 2 - function ip, then for 
sufficiently large n there exists j(k) such that ||0 n ||fc < npw 

Furthermore, one has some estimates for the first eigenvalue for closed mani- 
folds of dimension m > 2. Namely, if M is such manifold and its Ricci tensor 
satisfies 

p{X,X) > (m- l)a > (68) 
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for all X e TM and some positive real constant a, then for the first eigen- 
value of the Hodge laplacian we have (by Lichnerovicz formula) Ai > ma. 
On the other hand, for non-negative Ricci curvature, there is a lower bound 

7T 2 

Ai > -rp, where du is the diameter of M. With this newer result, the above 

M 

estimate could be improved as follows: 
Let M be as before and suppose 

p(X, X)>{m- l)a, aeR (69) 

then there are two cases 

(m — l)a it 2 
*i>- — T^ + IT */">o 

4 d M 

(m — l) 2 a dim) „ , . 

Ai < ^ — + ifa<0 (70) 



2 

M 



In the latter (3{m) is a real constant, depending only on the dimension. 
In the special case m = 2 we also have restrictions concerning multiplicities 
fij of the non-zero eigenvalues. For example, we have fij < 2 j + 1 for the 
sphere, fij < 2j + 4 for the torus, for the Klein bottle and the real projective 
plane the upper bound is 2j + 3, and finally fij < 2j — 2\{M) + 3 for surfaces 
of negative curvatur^§. 

Note that in the case of a Riemannian surface of genus g we also have the 
classical estimate for eigenvalue multiplicities: 

dim£ Afe = ±(2g + k+l)(2g + k + 2) (71) 



here and below x(M) denotes the Euler number of M 
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5.2 The ^/-calculus and Quantization 



5.2.1 The Fourier Transform, ^DO's and Sobolev Spaces 

The Fourier transform T is a well-defined endomorphism of the space of 
Schwartz functional S(R n ): 

v/ e s(R n ) f( P ) = Hf) P = [ /(■'■)< ,:i - r: <i"-r 

and with its inverse 

F- l {<j>) x = (27r)-» / (f>(p)e^d n p 

JM. n 

it is an involution for even functions due to the property T 2 {$) x = 0(— x). 
The fixed point of this involution is easily seen to be the Gaussian distribution 

f(x) ~ (27r)-"/ 2 e-l^ 2 / 2 

as one may easily check with a simple calculation. 

The map (27r)~ n / 2 .F is a unitary automorphism of the Schwartz space of 
function due to the Perseval-Plancherel theorem 

\\f\\u = {^) n/2 \\f\\u (72) 

No matter that S is the maximal closed under JF subspace of L 2 , it can 
naturally be prolonged over its dual space - the distributions in S"(lR n ) by 
the formula 44 ! 

{Fu,v) := (UiFv) (73) 
For example T(8) = 1 for the Dirac 5-function etc. 

One of the main advantages of the Fourier transform for the theory of PDE 
is that it turns differentiation into multiplication. More precisely, for every 
multi index a we have 

HD a f) =Cf D = i A (74) 

43 these are by definition the smooth functions / for which there exist constants Cki such 
that for every choice of the multi-indices, the estimate \x D f \ < Cm hold for all k, which 
means that the Schwartz functions, together with all their derivatives, decrease faster at 
infinity than the inverse of any polynomial \D l f\ < Cki(l + M)~ fc 

44 unlike ( , ) that we use for the usual dot product in R", here ( , ) denotes the scalar 
product in L 2 
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which allows for expressing a constant-coefficient differential operator by 
means of a polynomial in the ^-variable: 

P = J2a i D i = Op (P) :=T- l PT (75) 

Recall that the polynomial cr(P) = P = ^2 di is called the symbol of P 
(often denoted P (x, £) for simplicity) and we write Op(P) = P (x, D) for an 
operator with varying coefficients whose symbol is P. 

The notion of symbol allows for expressing easily the inverse of various dif- 
ferential operators. For example, the resolvent of the constant-coefficient 
operator 1 — A is simply given by 

For the ^-dependent case there is a natural extension of this formalism as we 
shall see later. 

Another familiar example is the resolvent of the heat equation, which is of 
particular interest for our considerations. As it not hard to prove by a simple 
Fourier technique, it has the form Op (e*'^ ), or as we shall formally write, 
e" tA . 

One may alternatively define the action of an operator P, acting on S(R n ) 
by its integral kernel K p (x, y), given by the symbol via Fourier transform 

P (/)»=/ K p (x,y)f(y)dy=(2it)- n [ e«*-Mp(x,t)f{y) dyd^ (76) 

We remark that the integral on the right side does not necessarily converge 
absolutely, so in the generic case we are not allowed to interchange the order 
of integration. 

Note that for integral kernels the sensible composition law is given by the 
convolution product 

(f*9)(x):= J f(x-y)g(y)dy= J f(x)g(y - x)dy, f,geS(R n ) (77) 

which is a kind of generalization of the usual matrix multiplication for con- 
tinuous indices. 

45 sometimes called also star product 
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Although we define the star product for the Schwartz class it extends to a 
broader functional space. For example in S , (M n ) the convolution algebra does 
not posses unit element. However, as we know, the ^-distribution is such an 
element in the extended case. There are some basic properties of * we shall 
take advantage of in the following: 

S*f = 

f * 9 = f-9, f-9 = f* 9 

* A = A * 

f*(g*h) = {f * g) * h (78) 

From now on we shall use both the symbol and the kernel of the operators 
in interest. 

As we saw in the examples, the inverse of a differential operator is no more 
differential, but belongs to a broader class of operators. For the study of 
these classes we investigate the spaces of symbols. Therefore we introduce 
the notion of the Sobolev space H s . For integer values of s, this is simply 
the space of functions with square integrable derivatives up to order s. The 
norm in this case is naturally expressed by means of the L 2 -norm with the 
formula 

ii/n.= £ip a /ik, 

\a\<s 

However, the definition is extended to all real values of s as follows: 

/ G H s ^/(l + |£| 2 )§ gL 2 

:= / (l + left'l/COPde (79) 

Actually, the first follows from the definition of the norm, thanks to the 
Perseval-Plancherel theorem. 

Of course, for integer s both definitions agree and of coursd^l H° = L 2 . 
Moreover, each two positive real numbers s < t we have the inclusion 

S(W n ) C H~* C H~ s CL 2 CH S CH ! C S'(M n ) (80) 

In this hierarchy of spaces we have quite a natural descending map - deriva- 
tion. Namely the derivative D a extends to a continuous map 



46 therefore from now on we shall write II . I In instead of 



o msueau 01 || ■ ||l 2 
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Intuitively, by applying this map, we loose the first \a\ square-integrable 
derivatives of / and fall down on a lower stair in the hierarchy. 
There are several nice properties, concerning the topology of the Sobolev 
spaces, which we omit for the sake of simplicity, but instead we give a good 
interpretation of H~ s as the dual of FT with respect to the L 2 scalar product, 
which extends to a perfect pairing between FT and H~ s . 
Baring these things in mind, we are ready to give the definition of a pseudo- 
differential operator. First let us consider a function p (x, £) : M n x W l — > C, 
such that 

1. p(x,£) is smooth in both arguments and compactly supported in the 
a;-space 

2. Va^GN" 3C a , p , \D°dIp(x,£)\ <c a ^i + \z\) d -w 

We call such a function a symbol of order d (p £ S d , d £ M) and assign to it 
an operator in the familiar way 

P(x,D)f = Op(p(x,0) = (^r n/2 I e« x -*Mp(x,t)f(y)dydt 

Such P we call the pseudo-differential operator {oi order d), associated with 
(or having as its symbol) p(x,£). 

First of all we need a few notes. The well-known differential operators are 
a private case of ^DO's according to the above definition - the ones with 
polynomial symbols. Their order as differential operators (the degree of the 
symbol) certainly coincides with the order given by our new definition. For 
a wider class of operators, however, the order is not bound to be an integer. 
For example, given a compactly supported smooth function f(x) £ C^°(lR n ), 
for all real d the symbol p(x,£) = f(x)(l + |£| 2 )2 £ S d defines correctly a 
^DO of order d. 

Another feature of the pseudo-differential operators is that they generalize 
the role of the PDO's in the sense that each ^/DO of order d extends to a 
continuous map 

P : H s -> W- d 

for all real values of s. 

And at the end of this section we give one more very important definition. 
We say that P(x, D) is smoothing if its symbol belongs to S d for all d, and 
that is C°°, but we write p £ S~°°. One can show that in this case P{u) is 
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smooth, regardless of weather u is smooth or not. 

5.2.2 Symbolic Calculus, Weyl Quantization and FIO's 

First of all we note that the operators under consideration constitute an al- 
gebra with respect to the operator composition and taking the adjoint. The 
infinitely smoothing operators form its ideal and in practice we mod out by 
it to obtain a nice factor-algebra. 

We use the notation ? 6 to show that p £ S d and = (~) d ty d will de- 
note the space of infinitely smoothing operators. We also use P ~ Q to show 
that operators (as well as the associated symbols) fall in the same equivalence 
class - differ only by infinitely smoothing term. We have the formulae 

a 
a 

whenever both sides are well-defined. In particular, if we consider only the 
leading symbols (the highest-order homogeneous parts of the symbol func- 
tions), we have simply 

o- L (PQ) = o- L (P)a L (Q) 

Moreover, the algebra thus defined, is graded, in the sense that if P G f i 
and Q £ \l/ m , then PQ £ ^i+ m - The above remains true also in the case of 
matrix-valued symbols. 

This construction may be used for example to construct asymptotic series of 
operators. More precisely, if pj £ S dj , dj — > — oo we have p ~ ^^Pj which 

3 

is unique modulo S~°°. 

Pseudo-differential operator possess the important property of pseudo-locality 
which means that if / £ H s is smooth in some open set U then Pf is smooth 
there as well. 

In practice one has to deal with fDO's defined on manifolds, rather then on 
M n as we implicitly assumed so far. The generalization, however, is straight- 
forward using partition of unity as long as one keeps in mind that the Sobolev 
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spaces the operators live in are preserved under diffeomorphisms and the 
leading symbols transform in a covariant way. This gives a sufficiently nice 
transition between coordinate charts, so that modulo ^-oo everything works 
fine. 

Now we mention some important classes of $£>0's that appear naturally in 
calculus and mathematical physics. 

Elliptic and Fredholm Operators Generally we call P G ^d elliptic if it 
is invertible modulo a ^-oo operator. More precisely, this would mean that 
there exists Q G ^ -d such that PQ ~ id and QP ~ id. Elliptic operators 
have some very nice properties compared to the generic case. One of them 
is the so-called Garding's inequality, which states that if P G is elliptic 
in a region U and / G Cq°(U), then the estimate 

\\f\U < C(||/||o+ H-P/llo) 
always holds for some constant C. 

Note that the sum of two elliptic fDO's need not be elliptic. This is always 
true, however, if they have both positive symbols. This allows for defining 
an elliptic ^/DO on a manifold again using a partition of unity. A crucial 
point in the theory of elliptic operators is the fact that when defined on 
closed manifolds, they are Fredholm. Fredholm operators are usually defined 
as bounded linear maps between Hilbert spaces Fred(L, M) G Hom(L, M) 
that are invertible modulo a compact operator, which means that if F G 
Fred(L, M) there exist operators G,H G Hom(M, L) such that FG — id is 
compact in M and HF — id - compact in L respectively. Such operators 
have finite-dimensional kernel and co-kernel (the null-space of the adjoint 
operator). Moreover, their image is closed (and the same refers to the image 
of the adjoint). It is not very hard to see that each operator with these 
properties is Fredhom and that they survive under composition and taking 
the adjoint. The above properties of Fredholm operators are widely used - 
for instance, for the derivation of Hodge decomposition in algebraic topology, 
or for the definition of index. 

k-$Z)0's For a natural generalization of the notion of a pseudo-differential 
operator, that is often used in physics, we first define the /c-Fourier transform 
by the formula 

W(0 = MO = [ e-^u(x) d n x (82) 
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and to each symbol function 47 ! p(x. <f). we associate an operator in the familiar 

way 

P[u] = Ov(p)[u]=T^\p]T h [u] 
In more explicit terms what we have is 

P[u]{x) = (±J ' J e*Mp{x,t)* k {t)<T£ ( 83 ) 

The additional fc-dependence appears in various problems as we shall be able 
to make sure very shortly. Apart from the /c-dependence hidden in the Fourier 
transform, one may have to deal with a more explicit one - in the symbol 
itself. In this case it is a standard situation to observe such dependence also 
in the estimate for p from the definition. Namely, one may be interested in 
symbol functions for which 

Va^GN" BC^, \D«DPp k (x,Z)\<C a ,pk r (l + \Z\) d -W 

and then one has an additional r-parameter which appears to be very helpful. 
For example, if the above is true for a whole series pj with corresponding rj 
than monotonically tend to minus infinity, then this series is an asymptotic 
expansion of some symbol pk in the sense that 

The above has the meaning of a power series expansion where the reminder 
is bounded and the bound is governed by a term of order k~°°. 
One significant subclass here is the class of the so-called classical symbols for 
which the above expansion is in integer power j 4 ^! of k: 

3 

In particular, the leading symbol of the above is <7l = k r po(x, £). 

The space of operators, associated to such classical symbols is usually denoted 

47 in the sense, clarified above 

48 note that the absence of upper index on the right hand side means that the symbols 
Pj (x,£) fulfill only the standard estimate for a generic symbol function, whereas for p k (x, £) 
it is with an additional factor k r as previously introduced 
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by W d . We note that a lot of \I/.DO's appearing in mathematical physics are 
of this type. For such operators we still have a nice-looking composition 
rule - let Pi,P2 £ and Pi,P2 denote their symbols respectively. Then if 
Q = P1P2, we have for its symbol q = u{Q) the power series expansion 



q{y,v) = exp 
evaluated at x = y, £ = r\. 



k 



(D x ,Ds) 



Pi(y,£)p2(x,v) 



(84) 



Weyl Quantization We already saw that (at least formally) a ^>DO is 
defined by its symbol (and vise versa) via the formulaF^l 



P[u](x) 



k_ 
2^ 



y 



In some cases, however, it is more convenient to use a slightly different defi- 
nition of a symbol, we shall refer to as Weyl symbol and denote p w , in order 
to have the above in the form 



P[u](x) 



k 

2^ 



V 



x + y 



,n u(y)e lk{{x ~ y) ^d n id n y (85) 



Note that this does not introduce a different operator, but simply offers 
another formulation. The relation between the two symbols is given by the 
formulate 

"1 



P (x, = exp 



P(x,0 



(86) 



Note that an alternative way to define a ^>DO is by its integral kernel. 
Namely, if K p (x, y) a matrix-valued function with compact x-support and 
f(y) is vector- valued, with compact support, then we define the operator 



Pf(x) 



K p (x,y)f(y)dy 



49 we refer to the more generic case of a k — ^DO, the usual one is, of course, the same 
with k = 1 

50 here and below we use the notation D = — iV 
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as a ^DO that is \IZ__oo if K p is smooth. For the case of Weyl symbols one 
may use the Fourier inversion formula to find the one-to-one correspondence 
between the two representations. It is obvious from the above that 



and then, by Fourier inversion 

One of the main advantages of the Weyl symbol is that (in terms of operator 
multiplication) it is written in explicitly symplectic-invariant form. In order 
to see this one needs to regard first (x,£) as a point z in the cotangent bundle 
over the x-space, endowed with a natural canonical symplectic form oj. Then 
if PQ = R, for the corresponding symbols it is true that 

r(z) = p w (z - \u>D„)q(z) = q W (z + \uD z )p{z) 

which can be interpreted in the following way: let 7 G Sp(n, R) and 7* denote 
its natural pull-back to functions 7*/ (z) = f('y~ 1 z). It is then a matter of 
computation to show that 

7*p (x, D) 7*g (x, D) = 7V (x, D) 

Using this invariance we can introduce a quantization scheme that is mostly 
used for fDO's. Let OY>(pf) G ^% j = 1,2 and let Q = PiP 2 denote their 
product with Weyl symbol q w . Then the latter is still a classical symbol and 
its asymptotic expansion is given by the formula 

00 1 / i \ k 

fc=0 ' ^ ' 

where the symbols p^ k and p^ k are those that take part in the asymptotic 
power series representing pY and p^ respectively, and {.,.} denotes the 
Poisson bracket. 
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One may readily see that the above formula leads to a natural correspondence 
(modulo terms of order k~ 2 ) between the classical quantities (functions on 
T*X) and their dynamics, defined by the Poisson bracket on the one hand, 
and the quantum observables (operators, acting in Hilbert space) with the 
commutator algebra on the other. More precisely, from the above we easily 
derive the formula 

a[P, Q] = ^W(P), <r(Q)} + 0(k- 2 ) (88) 

The Weyl quantization is quite convenient physical considerations as it de- 
fines the operator by its natural action on plane waves. Moreover, it allows 
for an easy expression of the expectation value 

(V | Op( P w ) I V) = (J^J J^p w (x,OW^d n xd^ 

where the Wigner function is defined as 

W^x, = / ${x - y/2)*p(x + y/2)e-^d n y 

A Glimplse On Fourier Integral Operators We note that the set of 
all k — fDO's forms an infinite-dimensional Lie algebra with the grading 

and in particular ^I/ 1 forms a sub-algebra. It is interesting to see what are 
the corresponding Lie groups (if there are any) to these algebras. It turns 
out that these groups exist and are formed by invertible /c-Fourier integral 
operators. To introduce k — FWs one takes the standard construction of a 
classical k — ^DO via an oscillatory integral and extends a little. Namely, 
the action of a k — FIO is given by an expression of the form 

Fk[u](x) = (-p\ [ e lk ^y>Vf(x,Ou(y)dydt (89) 

where / G S r d is a classical symbol in the above sense and the non-degenerate 
phase function <p is homogeneous map of degree one in £ and generates a conic 
Lagrangian manifold A e T*M x T*M/{0}. We may also restrict ourselves 
to the specific case where A is a graph of symplectic diffeomorphism 77 and 
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- its generating function. With this interpretation in mind it is natural to 
write the above oscillatory integral in the form 



With the choice S(x,^) = (x,£) one recovers the definition of a k — tyDO. 
In this case, as it is not hard to see, i] = id and A is the diagonal in T*M x 
T*M/{0}. Hence, we have the imbedding 



Now let us just enumerate some of the nice properties of the FJO's that 
make them an often preferred tool for calculations. 

First of all, due to their invariance under diffeomorphisms, these FIO's can 
be defined on manifolds^*] as bounded linear operators. Moreover, they ex- 
tend to such maps between Sobolev spaces, just as the fDO's do. Second, 
these operators are not pseudo-local, but posses the more general property 
to make the wavefront evolve in a fairly 'nice' way, namely 



It is also a important to see how these operators behave under composition 
and conjugation. The answer is quite satisfactory: 

1. Let Fi, i — 1,2 be two FIO's of order fcj and y), 4>i(^i,x,y) 
denote their amplitudes and phase functions respectively. Then F\ o F2 
is again a FIO of order k\ + k 2 with amplitude bi(C,2,x,y)b2(^,2,x,y) 
and phase function 0i(£i, x, y) + ^2(^2, %, y)- The corresponding diffeo- 
morphism is just a composition of the morphisms, defined by F\ and 
F 2} namely r] 12 = rji o rj 2 . 

2. Now let F' denote the conjugate of F. Then it defines a diffeomor- 
phism i]' 1 and is given by an amplitude b*(£,,x,y) and phase function 
—<f)(!;,x,y). From here we see that the product of a FIO with its 
conjugate is actually a ^/DO. 

51 in order to avoid confusions, let us restrict our attention to compact ones for the time 
being 




(90) 



DO C ^DO C FIO 



WF(F[u\) C Ao 



WF{u) 
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We give clear account of the fact that the above exposition on pseudo- 
differential calculus and applications hardly covers the minimum for any 
decent survey. However the hope is most of all to draw the reader's at- 
tention to the topic and to the matter and ensure a sufficient platform for 
further exploration. Our basic guides in the second half of this section are 
[55] as well as the appendix of the dissertation [31J. Both works have the 
advantage of clear and plain exposition. A more thorough investigation could 
be found in Hormander's monograph on partial differential equations. 
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6 Trace Formulae 



Throughout the present work we became familiar with some basic concepts 
in modern (and not that modern) mathematical physics, starting with inte- 
grability, quantization and arriving at the theory of PDS's and \l/-calculus 
at the end. No matter how significant these things are for the whole theory, 
we have to admit that nothing yet has been said about spectral geometry 
itself. 

The trace formulae introduced below are the basic tools of spectral geometry 
and reveal the essence of its philosophy. Only in this section we have system- 
atically exposed the theory of the correspondence between spectra, geometry 
and dynamics. 

To start with, we recall the classical Poisson summation formula: 

E /(*) = E / m e2m{p,k)d p = E /» ( 91 ) 



The above formula is crucial for the theory of theta and zeta functions, and 
it relates the sum over all lengths of closed geodesies (the right hand side) to 
the spectrum of the self-adjoint operator \fK on the torus TP (on the left). 
Thus, it inspires a whole series of formulae, involving operator traces and 
geometric quantities to be discussed below. 

Even the simplest case that we consider here can be very helpful in practical 
computations. It is worth nothing for example to show the extremely valuable 
identities concerning traces on the unit circle 



E^ - \hH 



x 2 — k 2 x 
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= — cot TlX 



^5(a; — k) = 1 + cos 2nmx 



m=l 



One of the first generalizations was made by Selberg, who obtained a trace 
formula for a particle, moving freely on a negatively curved surface X. The 
formula relates the spectrum of the Laplace-Beltrami operator over such 
hyperbolic surface to the length spectrum of the closed geodesies over it. 
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For a smooth, even test function h on X, with compactly supported Fourier 
transform h, the formula reads 



oo 




2 sinh (H 7 /2) 



(92) 



3=0 



where if* is the set of primitive closed geodesies and £ 7 - their lengths. 
The quantities pj are related to the spectrum of the Lapacian by pj = 
\j + k/4, Xj G Spec(A), but we may take the scalar curvature with the 
normalization k — — 1. Then the left hand side is, according to Lichnerovicz 
formula, a symmetric function of the Dirac operator with eigenvalues pj. 
Selberg trace formula appears to be a powerful tool for relating classically 
chaotic motion to quantum spectra and apart form that has proved useful 
in scattering theoryr 2 I.One may attempt for example to restore the whole 
geometry of a hyperbolic surface X from the spectrum {Xj}- The idea is 
to determine one by one the lengths of the primitive closed geodesies on X, 
which are known to be isolated and represent the fundamental cycles gener- 
ating 71"! (X). Once 7i"i is determined, X is easily represented as a quotient of 
the Poincare unit disk. 

There is one rather algebraic formulation of the above trace formula, that 
relies on the notion of induced representations. Let G be a Lie group, 
T - its subgroup and let L : T — > H(L) be a representation of T. We 
use F to denote the space of functions / : G — > H(L) with the property 
V7 G T, f(j.g) = L(j)f(g). Moreover, the space F is equipped with a scalar 
product induced by the scalar product in H(L) 



where p is the Haar measure on G. 

Now the representation we call induced by L and denote U L , appears natu- 
rally as 



52 recall that the monodromy matrix constitutes a representation of SU(1, 1), which us 
a principle bundle over the Poincare disk 





U L @U 



M 
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From here it follows immediately that if L is reducible, so is U L and vice 
versa. 

Now let U L be completely reducible U L = @jUjMK Then we have for its 
trace 

TrU L = ^^TrM j = 

3 

The Selberg trace formula then concerns the trace of the induced represen- 
tation. It states 

TrU L = [ [ f{g. 1 .g- 1 ) XL { 1 )d 1 du{g) (93) 
Jg/t Jr 

Here u(g) is the normalized measure on G/Y induced by \i. 

In the trivial case when Y = {e} and G is abelian, we recover the classical 

Poisson summation formula 

/ m di=Y, fix) 

6.1 The Heat Kernel and (^-regularization of the Effec- 
tive Action 

One of the best known examples of trace formulae in mathematics and physics 
is the so called Weyl expansior^ for the asymptotic fundamental solution of 
the heat equation. Here we are going to briefly introduce this technique and 
give some examples of its applicability from quantum gravity. 



/ M{f)dn(M) 
Jg 



6.1.1 High-temperature Expansion 

We consider the heat equation 

(A + |), = 

mainly for closed manifolds but also as a boundary-value problem. In both 
cases one may construct a parametrix, or a fundamental solution modulo 
^-oo, given by a power series in y/t with locally computable coefficients 

53 sometimes called also heat content or Seeley-DeWitt expansion 
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a,k(x), depending only on the curvature and its derivatives. 

To start with, the kernel of the resolvent for R n is known to be of the form 

K (x,y,t) = (47rt)~t e^^r 1 - 

For a curved manifold, however, the Laplacian undergoes certain deforma- 
tion and so does its resolvent. Nevertheless, some basic features remain 
unchanged: namely, the homogeneity in \x — y\/Vi, the translation invari- 
ance and the exponential damping for large 'distances'. 
Now let us consider a perturbative expansion of the resolvent kernel for a 
generic closed manifold. It is quite obvious that the Euclidian laplacian is 
the zeroth approximation in 

A = A + eAi + e 2 A 2 + . . . 

and respectively K° is the zeroth order term in the solution. Substitution in 
the non-perturbed equation yields 

(A + ^)K = 6 + R (94) 

where R is a lower-order kernel. We claim that the first correction ternj^l is 
given by Kq — Kq * R. This is easy to verify by a straightforward substitution 
in the initial equation - here we take advantage of the fact that * commutes 
with the heat operator and 5 * R = R by the very definition of 5. Thus we 
have cancelation of the first-order residue R, which is replaced by a lower 
order term —R * R. To get rid of it, following the same idea, we correct the 
Green function by adding the term K2 = Kq * R* R, which generates on its 
side another residual term of lower order. With this iterative procedure one 
obtains a convergent series of the form 

00 

K ~ J2 K i> K 3 := (-l) j Ko ?R---*Il (95) 

fc=0 j 

and as j tends to infinity, the order of the residual term in the approxima- 
tion tends to minus infinity, so the above series represents by definition the 

54 here * stands for the convolution product for the integral kernels defined in the usual 
way 
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parametrix of the heat equation for a generic smooth compact manifold with- 
out boundary. We note that this construction is in fact quite natural - this is 
nothing but a convergent geometric series in Banach space, where the usual 
multiplication is replaced by convolution of Voltera kernels and the condition 
\R\ < 1 makes sense once the corresponding operator norm is introduced. 
Let us consider now the kernel K(x, y, t). As the general theory says, it is a 
well defined smooth function everywhere in its domain, except on the diago- 
nal x = y where it has a singularity at t = 0. The idea of the Weyl expansion 
is to give an asymptotic representation of K(x,x,t), or more precisely of its 
parametrix, in the limit t — > as a power series. This is being done with 
the help of the theorem of Minakshisundaram and Plejiel, stating that the 
parametrix of K is really expandable in the form 

K(x, x, t) ~ V ~a k {x)t k - n ' 2 (96) 

(47T 1 ' 



2 

fc=0 



modulo infinitely smoothing terms with d being the dimension of the under- 
lying manifold. 

One immediate consequence of the MP theorem is that the coefficients in this 
series are spectral invariants, that locally depend polynomially, as we shall 
see later, only on the curvature and its derivatives. The first and simplest 
ones have very clear geometrical meaning!. In the case of closed manifolds 
for example, we have 

a = Vol(X) 

ai = - J Ry[gd n x 

a 2 = ^J x (2\R\ 2 -2\Rtc\ 2 + 5R 2 )^d n x (97) 

In the latter |i?| 2 and |-Ric| 2 are the square norms of the curvature and Ricci 
tensors respectively and R, as usual - the scalar curvature of M. Now we 
can easily see that the first approximation in high-temperature (or ultravio- 
let) limit gives us the opportunity to 'hear' the volume of the under lying 
manifold, as Mark Kac formulated it. The second is nothing but the Einstein- 
Hilbert action. In dimension two it is proportional (via Gauss-Bonnet for- 
mula) to the Euler number of the manifold. The third accounts for the 



55 



here and below we adopt the notation at = / a,k(%) for the integrated invariants 

Jx 
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presence of anomalies in AD field theory, as we shall explain below etc. 
Obtaining these quantities in explicit form however is not that easy. There 
are several preferable ways to do this - some are based on geometrical invari- 
ants and dimensional reduction, others use residual calculus for the resolvent 
kernel. Here we only give a brief idea of the most popular methods that have 
found application so far. 

Perturbative Ansatz To begin with, let us consider a Riemannian mani- 
fold with almost flat metric. One may regard the MP theorem as a construc- 
tive tool for perturbative calculation of heat invariants. Supposing that the 
curvature is relatively small, one may expand the curved metric as a sum of 
the flat (Euclidean or Lorentzian) metric rj^ and a correcting term 

Then one needs to obtain the variation of the Laplacian and by substituting 
the Euclidean heat kernel K in the formula (|95|) for the perturbed Laplacian, 
obtain (at least the leading part of) the remainder in (I95p . Next, using the 
convolution procedure described above, one may, in principle, find higher- 
order corrections as well. The coefficient in front of the j th power of t in the 
cut-off Voltera series so obtained is certainly a polynomial of the entries of 
h^ v and their derivatives. These polynomials may look confusingly complex 
(the whole complexity coming from the cumbersome-looking variation of the 
Laplacian) and yet they are not that difficult to transform in the form cij 
given above by means of Stokes theorem and some basic identities from vector 
calculus and differential geometry. The integral invariants are then obtained 
by simple integration. For example the coefficient ao is trivially given by 

a = (47rt)%TrK(x,x,t) lt _^ = Vol(M) 

However, analogous computation even for becomes complicated. Thus in 
order to obtain in relatively simple manner higher order coefficients, we need 
some more efficient method that are described below. 

Variational Method In this paragraph we shall be using the so called 
smeared heat kernel 

K(f,D,t)=Tr h2 (fe- tD ) 
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that emphasizes the distributional nature of K. Here / is a smooth function 
and D - an elliptic, self-adjoint operator, acting on the space of square- 
integrable functions L 2 . Accordingly, we refer to ak(f, D) as to the 'smeared 
coefficients 

Then, as this trace is still well defined, one may write the following variational 
equations for the corresponding heat coefficients a,k(f,D): 



d k (l,e- 2ef D) 
d k (l, D-eF) 
a n/2 ^{e-^F,e-^D) 

where the notation x is used for — 

de | e =o 
functions. 

To prove the first equation we use that 



= (n-2k)a k (f,D) 
= a k -x(F,D) 

= (98) 
and /, F are two arbitrary smooth 



Tr[exp(-e- 2e/ tD)] = Tr[2/iDexp(-tD)] = -2t— Tr[/ exp(-tD)} 

and expand both sides in power series in t. The second one could be easily 
checked in a similar way. In order to prove the third identity, we construct 
the operator D e j = e~ 2e ^(D — SF) for which, by the first equality, we have 
a n /2(l, D e} $) = 0. Then we vary the latter with respect to S and since the 
variations commute, we can easily apply the second equality to prove what 
we have in the third. 

The strategy we use to construct heat invariants, first used by Gilkey (35], is 
the following: first, according to the MP theorem, we construct a k as a linear 
combination of all possible independent invariants of the corresponding de- 
gree, one may build from the metric and its derivatives (and the derivatives 
of the potential term, if there is one). Then we use the above variational 
equations and some other generic considerations to eliminate the free pa- 
rameters. For example, if n = 2, there are only two invariants - E and the 
scalar curvature R. Note that we can always eliminate the derivatives of / 
via integration by parts, which is no longer true in the case of manifolds with 
boundary. Moreover, all half-integer coefficients equal zero in this case, since 
one cannot construct an odd dimensional invariant integral if boundary is 
absent - integration by parts makes them all vanish. 

Now, in order to make things a bit more complicated, we let all our con- 
siderations take place in a vector bundle V, so that we may have additional 
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indices (gauge, spin etc.) which we denote by Greek letters, or omit, if pos- 
sible. Then let us write once more the Laplacian under consideration, in the 
generic form fl56l) . which now would be 

D a/3 = _ g ik V a y /3 _ 

where V and E are respectively a (matrix-valued) covariant derivative and 
an endomorphism (as chosen in (I5T1) ). taking values in a vector bundle V over 
X. We may think of E as of a kind of effective potential, as far as physics is 
concerned. In the physical example of a free massless scalar field in curved 
space-time, it is given approximately by 

E°f> ~ -l([/"($))^ - "~ 2 R5 a P 
2 V y " 4(n-l) 

For the construction of the non- vanishing a,k(f,D) we consider several quite 
special cases. Assume, first of all, that the underlying manifold has a direct 
product structure M = Mi ® M 2 and the operator D is given respectively 
by D = Di ® 1 + 1 (g> D 2 . Since the bundle indices are also independent, one 
may write symbolically exp(— tD) = exp(-tDi) ® exp(— tD 2 ) and use that 
the smearing function factorizes as well f(x\,X2) = fi{xi)f2{x2) to obtain 
finally 

a k (x,D)= ^ a p (x 1 , a q (x 2 , D 2 ) (99) 

p+q=k 

The above equation may be used to show that the dependence of the heat 
coefficients on the dimensionality d is fairly trivial - they are all simply 
multiplied by an overall factor (47r)~ n//2 . In order to prove this, one takes 
Mi = S 1 , for which we know, from Poisson summation formula, that the 
high-temperature expansion has only one non-zero term. More precisely 



and applying the product formula for a k one obtains the recursion relation 
needed. 



Therefore, it is sensible to write 



a 



(f,D) = (4vr)-"/ 2 J d n x^gtr v {fa } (100) 
ai(/,D) = (4vr)-"/ 2 J d n x^tx v {f{a x R + a 2 E)} (101) 
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and so on. We already know that a® = 1, and can prove it once more by 
using our new approach. The coefficient a<i = 1 is easily determined by the 
second variational equation. For a\, however, we need to include higher-order 
coefficients in the system of equation which would give in the end a\ = 1/6. 
When boundary is introduced, we have in general 




where h is obviously the induced metric and denotes the j th normal 

derivative at the boundary. The dimensional dependence of the constant 

coefficients is shown in complete analogy to the previous case. 

In the following we need asymptotic expansion on the interval x G [0, 7r], 

namely 

jr(i,-flg,^)~yj±i (102) 

where B ± denotes the Neumann (resp. Dirichlet) boundary conditions. For 
the modified Neumann (Robin) case we have 

K(l, -dl B S +) ~\[l t + \+S\[l+ V(S 2 ) (103) 

Since a is absolutely the same as before we omit it and compute the next 
two in the expansion. We have 

a 1/2 (/, = (47r)-V f d^xVhpftivf 

JdM 

a^D,^) = (4ny n / 2 ±y^d n ^tr v (fR + QfE) 

+ / d n ~ 1 xVhtr v {l3tfL kk + (3fd v f + ftfS) 

JdM 

where denote the components of the extrinsic curvature of the boundary, 
or the second fundamental form, defined as L(Y, Z) = (VyZ, d n ) (d n denotes, 
as usual, the unit inward normal to the boundary). The unknown coefficients 
(3 can be determined with the help of the above result for the interval. For 
example, it immediately giveS = — /3f = 1/4 and p± = 12 (/3J is 

56 when M is an interval, the boundary integral reduces to a sum of contributions form 
its endpoints 
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certainly meaningless). 

In order to determine (3f and (5f , we use the conformal variational equation 
(the first formula) in fl98l) . 

The variation of the second fundamental form L is given by the variation of 
the metric and has the form 

^- L kk = -fL kk - (n - l)d v f 
de\e=o 

Now, in order to keep the Robin boundary conditions conformally invari- 
ant (in the Dirichlet case this is always so), we demand that the boundary 
function S fulfills 

j- S = -fS + \{n-2)d v f 
de\e=o 2 

and after finding a bunch of explicit expressions for the other terms in (I104p , 
we arrive at a linear system for the /3's, satisfied by 



fi+ = f3~ = 2, (3t = -(3., 



.3 



For the sake of generality we suggest below an explicit expression for the first 
few heat kernel coefficients of Kd(1, D,t) 

a = / dimV" 



x 



'IT j 

ai/2 = -"V / dimV 
z Jax 

at = \ [ Tr{6E + RI v ) + \ I Tr(L u ) 
6 Jx 3 J dx 

a 3/2 = / Tr(96£ + IQR - 8Ric a a + 1L\L) - 

iyu Jax 

(104) 

For the computation of more complicated terms we refer to [33] and [35J. 
In [5] and [61] one may find alternative methods that allow in principle for 
computing infinitelyl^l many a k s. We also refer to [69] for a survey on the 
famous DeWitt procedure, that is not mentioned here. 



57 this should not be taken literally, as it would mean 'with infinite complexity', since 
both procedures are iterative 
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One more use of the heat kernel expansion can be found in obtaining semi- 
classical distribution functions of various systems from quantum statistical 
mechanics. For examples 

1 



'S 1 



1 M + C + (105) 



Vg3 



4tt/3 V 3 15 

1 'i + p + £ + ... 



(4vr/5) 3 / 2 V 6 

are the high-temperature asymptotical for the partition functions of the 
rotators with two, one and zero symmetric axes respectively. Using this ex- 
pansion we are able to compute the thermodynamical functions in Quantum 
Statistics and QFT in the high-temperature (respectively ultraviolet) limit. 

For example the inner energy is given by U = — — — — , the Helmholtz free 

dp 

energy - simply by F = — — , the specific heat - by 

Cv = m : r and so on - 

The quantities find also one very curious application - in the theory of 
integrable systems. As we already saw, each such system is equivalent to 
an iso-spectral deformation of its Lax operator. In some cases we can re- 
gard the Lax operator as a Laplacian acting on a curved manifold. Then 
the spectral invariants in the heat-kernel expansion for this operator are first 
integrals for the corresponding hierarchy of dynamical systems. In the case 
for the KdV equation for example, the Lax operator is known to be of the 
d 2 

form L = — — — - + u(x), which is the Hamiltonian of a one- dimensional scat- 
ox 2 

terer with potential u(x). The spectral invariants in this case generate the 
Korteweg - de Vries hierarchy via the first-order PDE system 

du (2n)\ d , r . . 

m = «°"<I"W < 106 > 

where [u] n denotes the n-jet of u over the point x. For n = we obtain the 
trivial u = const, for n = 1 we have u t = u x , and n = 2 gives the well-known 



58 recall that (3 ~ — , where T is the absolute temperature and k - the Boltzman constant 
kT 
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KdV equation u t = u xxx + 6uu x . 

Using the above formula, one may in principle generate non-linear PDE's 
from the same hierarchy - as many as heat kernel coefficients at are available 
and on the other hand, reconstruct these coefficients from the explicit form of 
the integrable equations in given KdV hierarchy. This interesting observation 
is left to the reader for further development. 



6.1.2 The (-function Approach 

In the following we suggest an interesting illustration of the way the MP 
theorem may lead to rather unexpected benefits for modern physical appli- 
cations. The particular case, considered here is still rather simplified, but the 
technique used, is very powerful and somewhat conceptual, with far reaching 
consequences, some of which we discuss here as well. 

In order to get started, let us first introduce the notion of spectral (-function 
of an elliptic self-adjoint operator. 

We remind that the definition of the classical Riemann (-function involves a 
Dirichlet L-series of the form 

oo 

c(s) = J2 n ~ s ( 107 ) 

n=l 

which is convergent for 5J(s) > 1. However, Euler managed to introduce an 
analytic continuation in the whole complex plane (excluding s — 1) with the 
following functional equality 

c(s)= r(f y J c(i-*) (los) 



For the proof we use the famous Poisson summation formula ( 1911) and then 
the Mellin transform for (: 

1 f°° 

((,) = -—/ t-^e-** (109) 

To obtain this formula one only needs to expand ((s)T(s), using their defi- 
nitions^ and then interchange the summation and integration. 



POO 

59 recall that T(s) = / i s_x e _t d* 
Jo 
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In this way (by means of formula fllOSp ) ( is defined as a meromorphic func- 
tion with simple poles at s = 1, 2 ... | for even n, and s = — ~, | . . . ~ for odd 
n, as well as simple zeroes for s G — N - the poles of the Euler T - function. 
There are also non-trivial zeroes, that are proved to be bounded in the strip 
< 3?(s) < 1 (none of them lies on the real axis) and spread symmetrically 
with respect to the line 9ft(s) = 1/2. The famous and still unproved Riemann 
hypothesis consists of the assumption that all of them lie on this line. 

One-loop Regularization Let us now take advantage of the nice ^-formalism 
developed above. Namely, in physics we use (^-functions for regularization of 
infinities. To give an example what this is, we consider once more the func- 
tional equality (11081) . In the case s = 2 we have the resuli@ 




and via (11081) this defines correctly the value 

C(-i) = -— 

which is quite an amazing result: we sum all natural numbers and obtain a 
negative fractional number. The key to understanding this is to accept the 
value —1/12 as the regularized value of the above sum, the value we obtain 
when we get rid of divergencies via analytic continuation. 
In Physics we resort to the trick of regularization quite often in order to 
compute measurable quantities - vacuum energies, functional determinants 
and self-interacting potentials are some of the things that often need to be 
regularized. In the above example, we obtained a regularized expression for 
the resultant force, acting on a particle in an infinite chain with potential of 
interaction U ~ — x 2 and lattice constant equal to one. 
In order to solve more serious problems, we need first to build a more powerful 
formalism. In complete analogy with the above case, we define the Spectral 
£ - function of an elliptic self-adjoint operator D with the formal series^! 

Cd(s)= Yl A " Sdim ^A (HO) 

X&Spec(D) 

60 again a result due to Euler - in fact one may relate Cflk to the non-zero Bernoulli 
numbers, using a consideration about the harmonic oscillator given later in this section 

61 later we shall skip the dim£\ - term, assuming that in the degenerate case multiplicities 
are always encountered 
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which is correctly defined at least for > n/2. The question of convergence 
in the elliptic case is ruled by considered above estimates for the spectrum 
of D. Note also that the Mellin transform fl 1 9 j) applies here as well with 
the only difference that the sum over natural numbers is replaced by a sum 
over the spectrum^! and £(s) is so defined as a meromorphic function with 
zeroes, determined by the poles of T(s). 

For a more detailed study of the spectral ^-function and its analytical prop- 
erties we refer to [3D]. 

Now we shall focus on a more 'physical' application of £-regularization. 
Namely, in the following example we regularize the determinant of an el- 
liptic, self-adjoint operator, acting on functions in L 2 . The determinant in 
question appears naturally in the path integral formalism of - when we cal- 
culate Gaussian integrals for the partition functions in QFT. In our example 
we regard a free massless scalar field 4>(x) in curved space-time (X,g). The 
classical action is well known to be given in this case by the expression 

S = \J x C = lJ x ^^^VW\d n x (111) 



Variational principle yields V 2 = inside the domain of interest! 63 !. The 
kernel of the propagator is given by the path integral 

rif fiS\ 
K(q in ,q f ;t in ,t f ) = D p D q exp{—j (112) 

•J Qin \ / 

where q in = q(t in ) and qj = q(tf) are the initial and final positions respec- 
tively, 

dp k dq k 



D„ D p = lim 



n 



2nh 

k=l 



dp 
2nh 



- the integral measure, defined as usual, and S = I C(q, q)dt is the classical 

J tin 

action along the trajectory starting from qi n at time £j n and arriving at qf 
at time tf. It is clear, however, that apart from the classical one, there are 
infinitely many possible trajectories and K(q in ,qf;t in ,tf) is interpreted as a 
partition function in the space of such paths. 

In QFT due to the hyperbolic signature of the metric, the integral measure 



62 in other words, £(s) is the Mellin transform of the heat kernel 

63 and with the corresponding boundary conditions as long as boundary is concerned 
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above is somewhat ill-defined. Therefore the simple trick of analytic contin- 
uation in the complex plane is performed - the so called Wick rotation. This 
is nothing but rotation of the time coordinate from the real to the imagi- 
nary axis (t — ► it), so that we end up with an Euclidian or Riemannian, 
rather than pseudo-Riemannian metric and imaginary time. Due to the co- 
variant transformation law we obtain the transformed action (we refer to as 
Euclidean action from now on) as Se = iS and the path integral from above 
in the fornix 

Z[</>] = J exp (S E /h) = det D~ 1/2 (113) 

Let us write down the above partition function in the form e~ w . The func- 
tional 

W = ^lndetD = -^C'(O) 

in this expression we call from now on the one loop effective action for the 
corresponding theory. The effective action gives the one-loop quantum cor- 
rection to the classical one and is a powerful tool for various calculations. In 
our somewhat simplified considerations the path integral representing K is 
Gaussian (the classical action contains only a quadratic term) and we obtain 
immediately K ~ det(D)" 1 ^ 2 = exp(|C^(0)), so at least here W recovers 
one-loop contribution to the path integral exactly. In more complicated set- 
tings, where there are several quantum corrections to the effective action. 
Namely, we expand the Lagrangian up to the quadratic order in quantum 
fluctuations (f> 

L~L c/ + (0,J) + (0, J D0) (114) 

where J is an external current and L c i - the classical Lagrangian, to obtain 
the more generic form of the path integral 

Z[J] = e" Lci det (D)~* exp [^JD^J] (115) 

and in this case the effective action gives only the determinant contribution 
to the overall propagator. 

We are now close to obtaining a regularized expression for the effective action 
in terms of the spectral ^-function. What we have to do is simply express 
Cd{s) as the Mellin transform of the heat kernel and then differentiate at s — 



6 the measure is chosen in such a way that the final result is free of multiplicative 
constants 
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0. When we deal with ultraviolet divergences it is useful and legal to expand 
the heat kernel as a power series in t and neglect the higher order terms, as 
the MP theorem teaches us. Then, after introducing a cutoff r —>■ + and 
exploiting the formulae for the first spectral invariants, differentiation term 
by term gives: 



d A xy\g\ 
32vr 2 



+ — + ( — R 2 + — RnuR^ } In I -I 
2r 2 6r V 120 60 " 1 1 



where the remainder consists of finite terms only. Next we introduce a mod- 
ified action for the background field (gravitation) so that its singularities 
cancel the singularities of W: 



& gravid nv] ~ I ^X\/\g 



grav ij 



R + 2A b ( R 2 R UV R^ 

+ tt 6 — + 
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The quantities with subscript b are called bare coupling constants of the 
theory. They are functions of r chosen to cancel the divergent terms in W. 
Their presence can never be observed since the back-reaction can never be 
switched off. However, with their help, one obtains the expression of the 
regularized effective action as 

w reg = w + s b grav 

For the sake of honesty we note that there are several different ways to get 
rid of divergencies in physical theories, but we describe only the one that 
corresponds most directly to the topic of this work and at the same time is 
usually referred to as the most elegant regularization scheme known so far. 
We refer to [55] for a systematic exposition of the (^-regularization procedure 
from a physicist's point of view. Another good source of information for the 
application of this technique in QFT is [69]. The latter is in fact suitable 
reference for the most part of this section. 



The Oscillator There is a nice example that allows for exact expressions 
for the heat kernel and the operator (^-function. This example is one of the 
very few exactly solvable problems in quantum mechanics, but also one of 
utmost significance for the whole theory. 

The partition function of the harmonic oscillator can be written in the form 

1 1 °° R 
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with B 2 j being the non-zero Bernoulli numbers. 

For the relevant (-f unction, we have already a number of ways to retrieve the 
result 

Cho(s) = ^Us) 

with £r standing for the classical Riemann (^-function. 

One may find exact values of thermodynamical quantities in this and just 

a few other cases. However, in quantum statistics, the oscillatory degrees 

of freedom play central role, so the result, no matter how simple, is not to 

underestimate. 

Residual Calculus Here we give one more idea how to calculate heat 
invariants. The trick is to use the integral representation ( 11091) of Cd( s ) an d 
its inverse 

TiK(x,x,t) = — <p t-T(s)( D (s)ds (117) 
2m J 

where the contour of integration encircles all poles of the integrand. Then, 

d 

baring in mind that a k = TiK(x, x, t), we easily obtain the expression 

k\ ox k 

a k = Res s=n/2 - k T(s)( D (s) (118) 

and in particular a n / 2 = C(0)- Curiously enough, this gives £(0) = for 
odd-dimensional closed manifolds. 

The power-series form of the heat kernel, obtained above, allows for express- 
ing the high-temperature asymptotic and all thermal functions once we know 
explicitly the (^-function in certain quantum-mechanical problem. The latter 
is, unfortunately, very rarely so. At least, the above expression gives one 
more curious connection between the analytic properties of a Laplace-type 
operator, and geometric quantities of the underlying manifold. For example 
it is far non-trivial that the regularized range of the Laplace-Beltrami oper- 
ator on two-dimensional manifolds is proportional to their Euler number. 



6.1.3 Conformal Anomaly and the Polyakov Action 

Now we briefly discuss the case of conformal anomaly. In the classical con- 
formal field theory we have action, which is invariant under conformal trans- 
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formations of space-time, that are, roughly speaking, all isometries plus di- 
latation^] of the type 

Such classical theories are characterized by traceless energy- momentum ten- 
sor 

2 8W 

T„ = —P- (119) 

since the variation of the action with respect to dilatations is given by 

SW=- I 5pT^du g (120) 

However the conformal symmetry is broken by quantization even in the one- 
loop approximation^!, as one may easily verify, substituting the effective 
action, constructed above and obtaining finally 

T£ = a n/2 {x) (121) 

So, the one-loop quantum theory possesses an energy-momentum tensor, that 
is not traceless anymore. 

There is a very famous physical example in which the explicit form of the ef- 
fective action may be obtained via integration of the conformal anomaly. The 
model considers a free scalar field on compact Riemannian surface, described 
by the classical action 



S= / d 2 x^/gd u (j)d u <j) (122) 
Jm 

It is well known on the other hand that all two-dimensional metrics are 
conformally flat, so we may choose a gauge in which 

and calculate in this gauge the trace of the stress-energy tensor and thus the 
conformal anomaly. Varying S we obtain 

5W =^~ [ d 2 x5p^d^d vP (123) 

127T J M 



65 and more precisely speaking, a pseudo-Riemannan space with metric of signature (p, q) 
has SO(p+ 1, q+l) as a conformal group - for example the conformal group of the Gaussian 
sphere CP 1 is easily shown to be PSX(2,C) - SO(3, 1) 

66 recall that in this case the laplacian appears as a fluctuation operator of the theory 
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Baring in mind that in the so chosen gauge 



we obtain after integration the conformal part of the effective action in the 
form known as the Polyakov action in string theory 



The inverse of the laplacian is naturally obtained via integration over the 
Euclidean Green kernel. 

We note that this results may in principle be generalized to higher dimensions 
but, unless M is conformally flat, the integrated anomaly represents only a 
part of the effective action. For more thorough investigation on this matter 
and further applications of spectral methods in quantum field theory we refer 
to [9], [57] and [70] . together with the references therein. 

6.1.4 A Note On Index Theory 

Now we mention a rather peculiar appearance of the spectral invariants from 

the heat kernel expansion in index theory. 

To begin with, define the index of a Fredholm operator as 



From this definition it is easy to see that the index is additive (for a com- 
position of operators) and ind(F*) = — ind(F), therefore the index of a self- 
adjoint operator is always zero. 

The index has a very strong topological meaning in the elliptic case - it 
depends only on the homotopy type of the leading symbol of the operator 
within the class of elliptic ^DO's of the given order. This observation has 
some deep consequences connecting analytic theory of elliptic operators on 
manifolds to algebraic geometry. 

A good basic example is obtained if we roll up the whole de Rham com- 
plex, distinguishing only between even and odd dimensional forms. Then 
the index of the Dirac type operator P := d + 5 : f^ ' — > defined as 
indp := dimker(P) — dim coker(P) = dimker(P) — dim im(P) 1 - gives ex- 
actly the Euler number of M. 




(124) 



ind-F = dimker(F) — dimcoker(F) 
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There are several similar observations that inspired the famous index theo- 
rem of Atiyah and Singer, relating the analytic index with the topological 
index, given by the A- genual of the corresponding bundle. Namely, for an 
elliptic operator, acting on section of this bundle, we have 

indP= / (A) (125) 

J M 

Suggesting that the integral on the right is integer-valued is already quite 
non-trivial. The approach of Gilkey to this matter uses the construction of 
the two operators P + = PP* and P~ = P*P. It is easy to see that their 
non-zero eigenvalues coincide, so when we consider the difference of the their 
heat-kernels, only the terms, corresponding to the zero modes survive, and 
these are the terms, independent of t in the series expansion. Thus, what we 
have in the end is 

mdP = a n/2 (P + )-a n/2 (p-) (126) 

One may consider once more the previous example - the rolled de Rham 
complex examine the heat kernels of dS and Sd, or more precisely, the coef- 
ficients di in both cases, to obtain the same result. Moreover, by the above 
construction and Gauss-Bonnet theorem, we have 

md(d) = x(M) = — I Rda 
2vr J M 

which is a special case of (11251) . 

6.1.5 The ^-function and the Berry Phase 

We conclude with a short discussion on //-functions. Consider first an op- 
erator P that is not positive definite. Then the ^-function is in general 
ill-defined. We define a new spectral function instead by the series 

V p{s)= Yl si S n MH~ S (127) 

fi&Spec(P)/{0} 

For this new quantity we have an integral representation again thanks to the 
Mellin transform 

_|_ 1 POO 

r) P (s)r(——)= t £ 2 i Tr[Pexp(-tP 2 )]rft (128) 
2 Jo 

67 the j4-eenus is defined as A = IT - c , z / 2 , where the c,'s are the Chern classes of 

o l iz stnh(Ci /2) ' L 

the bundle 
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In particular, for the ^-invariant we have 

f°° 1 

r)p(0) = / -^Tr[Pexp(-^P 2 )]dt (129) 

that is a regularized expression for the spectral asymmetry of P (recall that 
the (^-invariant Cz?(0) gives in a similar way a regularized expression for the 
range of D). 

Another useful applications of ^-functions is the computations of functional 
determinants for operators of the type under consideration. First of all we 
use that P 2 possesses a well-defined ^-determinant and thus 

det|P| =exp(-iCp 2 (0)) (130) 

Next we easily decompose Spec(P) / {0} to purely positive and purely nega- 
tive eigenvalues in order to obtain a formal expression for Cp(0)- The final 
result is 

det(P) = exp[--(C p2 (0) - VP (0))] det |P| (131) 

The argument of the oscillating exponent to the right is sometimes called 
Berry phase of the determinant and plays an important role in path integral 
formalism. 

The spectral invariants considered in this section possess also a topological 
meaning due to certain index theorems. They are of different sensibility and 
some of them are rather subtle for investigation. A good survey on this 
matter is [6lj 



^- + V(q) + ihd t )iP(q,t) = (132) 



6.2 Path Integrals and the Gutzwiller Trace Formula 

In this section we consider the non-stationary Schrodinger equation 

2m 

where the first two terms represent the quantized hamiltonian H(p,q), con- 
sisting of potential V(q) and kinetic energy 2— = -^-V 2 . It is easier to 
expand the solutions over the basis of solutions of the stationary equation 
- the eigenfunctions of the Hamiltonian (p^ : Hip^ = E^^. For bounded 
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systems this basis is known to be orthonormal and complete. The solution 
i/j(q,t) we are interested in is then expanded as 



ip(q,t) = ^c n e ^"Vn(<?) (133) 

n 

where the Fourier coefficients are given by c n = / (f) n (q)'ip(q,0) d m q. We 
can write it also as 

#7,*)= / K(q,q',t)i;(q',0)d m q (134) 



>E„t 
h 



where the integral kernel of the propagator K(q, q', t) = 4>n{q')4>n{q)e 
carries the whole dependence on q and t so it is itself a solution of the initial 
Schrodinger equation. Moreover, it is a fundamental solution^! due to the 
completeness of the basis. 

The Fourier transform of K in the time variable is a well-defined distribution 
depending on the dual variable - the energy. This is nothing but the resolvent 
of H and is usually introduced in a regularized way as follows: 

G(q,q',E + ie) = - / dte-^ E+i£) K(q,q' ,t) 
in Jo 

n n 

Once we have guaranteed the existence of the above integral, we may get rid 
of the imaginary shift by taking the limit e — > 0. In the following we shall 
be interested in the trace of the green function 

TtG(E) = [ G (q,q,E)d m q = J2T^-Tr (135) 
which is connected to the spectral density d{E) = \^5(E — E n ) with the 

n 

famous Sohotski formula 

ird(E) = -3 (TtG(E)) (136) 



Hhat is K(q, q', 0) = S(q - q') 
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6.2.1 Semi-classics for the Schodinger Equation 



For a relatively slowly varying potential 69 ! we are allowed to apply the Wentzel- 
Kramers-Brillouin (WKB) ansatz for the wave function: 

if>(q,t) = A(q,t)e% R ( q 'V (137) 

We note that this h — > asymptotic solution cannot be prolonged trivially 
for large times due to the formation of caustics at time t ~ H°. Therefore, we 
need to let it evolve with clear account of the singularities it passes through. 
Substitution into the initial equation yields (after separating the real and 
imaginary part and taking the limit h — ► 0) 

d tP + v{ P v) = o 

d t R + H(VR,q) = (138) 

The former is interpreted as a continuity equation for the flow p = A 2 (q, t), v 
— Vi2(g, i), referred to as the Madelung flow, while the latter is simply the 
classical Hamilton- Jacobi equation for R(q,t). 

Now it is easy to obtain the evolution of R with respect to this flow, replacing 
VR -> p 

R(q, t) =pq- H(p, q) = C(q, q, t) 

where H and L are respectively the classical Hamiltonian and Lagrangian of 
the system. Thus for a classical trajectory, starting at the point (q',p'), R 
evolves with time as 



R(q(t),t) = R(q',0)+ [ C(q,q,t')dt' 

Jo 



the integral to the right, being the Hamiltonian principal function. Denoting 
it by R(q, q', t) we obtain the classical initial and final momenta 

p = -V q >R(q,q',t) 

p = VR(q,q',t) (139) 

The conservation of matter demands covariant transformation of the volume 
form so we have 



A(g,*) = t/det(|^)A(g,0)=e- 




detl^ 



dq 



A(q,0) (140) 



69 the semi-classical approximation is valid only in the case when the potential does not 
change considerably within the de Broiglie wavelength 
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where the quantity v(q, q', t), called the Maslov index, counts the sign changes 
of the Jacobi determinant along the path, connecting q' with q. Strictly 
speaking, there are usually many classically allowed paths 7&, ending up in q 
at time t, and each such path has its own starting point, principal function 
and Maslov index, and gives its own contribution to the evolution of states. 
Thus we need to sum over all such trajectories: 



iR k ( q , q ' k ,t) I f dq',\ 



Mh,o) ( 141 ) 



Now we are ready to express the propagator, but first, to get rid of the 
singularities due to the point-like contribution of initial states in configuration 
space, we switch to momentum space, where the Dirac delta-distribution is 
transformed to a constant factor by simple Fourier transform. We have 

p(q,t)d m q = \c\ 2 d m p' 

and changing from p' to q in the righthand side, we easily obtain 

p(q,t) = \c\ 2 det(-d q >d q R(q,q',t)) 

For normalization reasons, the constant is easily shown to be c = (2irih)^ m ^ 2 
so, expressing the propagator as before, in terms of the evolved states, one 
obtains the famous expression for the semi-classical Van- Vleck propagator 



K(q, q', t) = (2mK)- m ' 2 £ e !£ ^ "t* y^det (-d q ,d q R(q, q>, t))\ (142) 



6.2.2 The Method of Stationary Phases and the Green Function 

From the above formula, it is not difficult to obtain the Green function as 
G(q, q', E) = ^G k with each G k being the Fourier image of the k th term 

k 

in the above sum, representing K(q, q', t). First of all, we need to know how 
to compute integrals of the type 

1=1 u(x)e lkf{x) d m x (143) 

in the limit k — > 00, where u(x) and f(x) are sufficiently smooth and ^s(f) > 
0. The idea for this method is to take advantage of the fact that the main 
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contribution to the integral comes from the set of points where the phase 
function is stationary - that is /' = 0, and ^s(f) = 0. 
Around a stationary point xq we have, up to the leading term 

u ( x ykf( x)d m x = elkf ( X0 ) ( ^Y 1 ' 2 <*o) 0(k~ m ' 2 - 1 ) 

\kj Vdet(-i/"(xo)) 

(144) 

With this formula in mind it is not hard to obtain a good approximation 
for some of the integrals, representing Gk{q,q',E). The stationary phase 
conditional is given in our case by d t Rj- + E = since z/ fc (g, q', t) is constant 
on the trajectory. This means that the time of the saddle point t c is the time 
it takes to a particle with energy E to move from q' to q following j^. Thus the 
stationary phase method is good only for relatively long trajectories. A direct 
application of the above formula (or rather to its direct multidimensional 
generalization) yields 

o i r(q, q \ e) « g^p: x (145) 

x [det {%R(q,</,t e ))]~* J\det(-d q ,dqR(q,q',t c ))\ 



where S(q, q ,t) = R(q, q \t c ) + Et c = / pdq is the classical action func- 

Jq> 

tional for the trajectory q' — > q. 

We may, on the other hand, separate the spatial derivatives of R in normal 
and tangential parts with respect to the flow and thus represent the latter 
determinant as 

det (-d qi d q ,R) = -i- d?R(q, q', t) det (-d^s) (146) 

which leads an integral expression for the long-time contribution to the trace 
in a form, suggesting once more a stationary-phase approximation. Here the 
stationary condition is 

VR(q, q', E) = V q S\ q=q , + V q ,S\ q=q/ = p - p' = (147) 

which means, together with the energy-preserving condition, that the orbits, 
giving the main contribution to the trace are not only periodic in the spatial 



711 



the condition, under which the approximation is legal 
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variable, but also in momentum, so they end up in the same point in phase 
space. 

Introducing the stability (or monodromy) matrix P 7 of an orbit, as the matrix 
relation between the initial and final variations of the phase space variables 
in direction, perpendicular to the flow0 

8(q,p) ± = P J 6(q',p') 1 (148) 

and using the definitions of the initial and final momenta and their variations 
in terms of the action derivatives, one may prove the matrix relation 

det(Id - P 7 ) = det (d q d q S(q, q, E)) det {ti$ q S)~ X (149) 

which allows for the expression of the long-time contribution in somewhat 
completed form. Namely, after all substitutions and integrating over the 
diagonal q = q' we obtain 

iS k (E) iirv k (E) 

TrG l r(q, q', E) « I / ^ \ \ ~ ° = (150) 



iU J *\\ J|det (I-P+) 



where summation is taken over all closed classically allowed paths through 

I dq\\ 

Now realizing that the integral to the right is simply the period (p — — = 



1\\ 

dt = Tfc and taking into account the revolutions of simple orbits, that also 



give contribution to the tracqlj we finally obtain the famous Gutzwiller trace 
formula 

Tr(G) = Tr(Go) + \ Y T e^-^ (151) 

^tT|det(I-4)|l 

The contribution from short-time trajectories on the other hand, for which 
the stationary phase method is not applicable, is approximately given by 

short f .. .1 t?\ 1 [°°t m \f j(^9^-V(,)t+E, 



im 



71 note that this is by definition the matrix of the linear Poincaree map for this orbit 
if I is the winding number of the orbit, then we simply have T — * IT, S —> IS, v — > Iv 
and P-y Pl f 
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where S (q, q', E) = ^2m{E -V)\ q — q'\ is the short-distance approximation 
of the classical action. Moreover the Hankel function of the first kind is 
defined as 

Now one simply takes the imaginary part of the trace and substitutes in the 
integral the asymptotic expansion for the Bessel function 

W-^Vi)(£f, N«i 

which is the real part of part of H^\x). 

Then, after integration over a ball of radius \p\ = ^j2m(E — V) in momentum 
space, one obtains the famous Weyl law for the asymptotic distribution of 
energy levels 

,. Ar/X , Vol{M)\i 

hm N(X) = -r — 

A-oo 2 d 7rfr(f + l) 

6.3 The Berry- Tabor Trace for Integrable Systems 

In the derivation of the previous result, we implicitly assumed that the closed 
orbits are isolated so that the stationary phase method for the Green function 
could be applied to each orbit separately. In practice, however, this is not 
always a reliable assumption. There are cases when the closed orbits come 
in continuous families filling whole areas of the phase space. One such case 
is a bound integrable system, with multiply periodic classical motion. In 
this case, as we know, phase-space trajectories fill invariant tori, given by 
action-angle variables. The idea for the derivation of the following trace 
formula is to use these variables in order to simplify the trace expansion in 
the integrable case. 

We start with the familiar expression for the semiclassical propagator 

K(q, q', t) = (2nih)- m/2 ^ e ^%^ -^ ,J\det (-d q >dqR(q,q',t))\ (153) 

k 

where the summation in our case is taken over all closed paths, and without 
writing additional subscripts, it is understood that the principal function 
R(q, q', t) is evaluated separately at each such path. 
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Next, we use the fact that the action variables Ik are constants of motions 
and the angles 9k evolve linearly with time according to 9 = V iH(I) = 
so that we have the classical action in the form 

S = 1.(9 -9')- H(I)t 

and thus for the determinant in the propagator we can write 



det(-d q ,d q R(q,q',t)) = (^) 



1 



t m det(du/dl) 



When we take the propagator trace, we use the multiple periodicity with 
respect to the angle variables and expand the integral as a sum of integrals 
over the lattice, defining the m-torus with generators / and coordinates 9: 

r /*27r 

/ d m qK(q,q,t) = / K (9 + 2nL, 9, t) d9 

In fact, the integer vector L accounts for the different topologies of closed 
paths 13. Under the assumption that no vector L can give paths in more 
than one m-torus simultaneously, all summations considered so far can be 
labeled by the corresponding Z m - coordinate L. For instance, the classical 
action over a closed L-path is given by 

S L = 2irI h (t).L-H(I L (t)).t 

The above expression is independent of 9 and so should be the semiclassical 
propagator in these coordinates. Hence, integration over the angle variables 
is simply reduced to a factor (27r) m . 

After substituting it in the initial expression for K, we need to extract the 
zero-length orbits with L = - they give the smooth term in the spectral 
density, just as in the Gutzwiller case. To the rest we apply stationary-phase 
method and after that evaluate the amplitudes by introducing curvilinear 
coordinates £j in the /-space. Let £ point in the VH direction, and the rest 
m — 1 £j's be local parameters on the energy contour H(Il) = E. Then we 
finally obtain 

n(E) n(E) I - ' r"H.W»-t/4)HM 

< i ^ > fie*")/* £j |L|<"-i>/2M/ L )|/M7I)l 
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in the two-dimensional case for example, we sum over toric knots 
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where k is the scalar curvature of the energy contour in the ^-coordinates, 
/3l is the excess of positive over negative eigenvalues in the matrix 

d 2 I 

and the vector-valued Maslov factor z/l encounters the number of caustics 
met by the orbit at each cycle of the corresponding torus. 
The above proof of the Gutzwiller formula follows [IS] and for the Berry- 
Tabor trace we refer to the original paper [To] . 

6.4 The Wave Trace 

The wave trace is defined as the distributional trace or as asymptotic series 
for the kernel of the wave resolvent 

U(t) =Tr exp(zr^A) = e irVX (155) 

XeSpecA 

The very crucial fact that makes the above a powerful tool in inverse prob- 
lems, known as Poisson relation, is that the singular support (wave front) 
of U(t) is contained in the length spectrum (the set of lengths of closed 
geodesies) both for closed manifolds and manifolds with boundary^. More- 
over, as to each length in this set always correspond at least two geodesies 
(say 7 and 7 _1 , where the latter is understood as 'time reverse'), we may, 
with no loss of generality, restrict ourselves to the even part of the trace 

U+(r) = Tr cosrv/A (156) 

What we have to deal with in practice is the singularity expansion of the 
above trace around closed geodesies. More precisely, as long as all closed 
geodesies are non-degenerate, one has (modulo \lZ_oo) 

U(r) = a + Yl a ^ T ) ( 157 ) 

l£Lsp(M,g) 

where the distribution is given by a power series expansion in r with 
coefficients that actually coincide at r = with the heat invariants. The 

an interesting observation in the case of real-analytic manifolds shows that we may 
have an equality rather than inclusion if we consider the analytic wave front (defined as 
the complement of the set where U (r) is real analytic) instead of the ordinary one 
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more complicated terms ot\ involve infinite series with logarithmic terms, but 
the principal part that measures the most severe singularity - that is the 
coefficient in front of it — I + ie)^ 1 , is given by the familiar formula 



where the sum is taken over all closed geodesies with length Z 7 = I and Ijf 
denotes the primitive length in each case. This formula works both for closed 
manifolds and domains with Neumann boundary conditions. In the Dirichlet 
case one has to modify it with a factor (— l) r , where r is the number of 
reflections of the given orbit (note that this is the well known from geometric 
optics phase inversion from a mirror). In the Robin case this factor is more 
complicated and will be discussed below. 

From the above it is already quite visible that the wave trace is a more 
general case of the Gutzwiller trace formula that we derived in connection to 
the Schrodinger equation. The information contained in the latter will often 
be sufficient for our considerations, but knowledge about the full wave trace 
(if possible at all) is much more profitable in inverse spectral problems. 
There are several distinct approaches aiming the computation of the wave 
trace coefficients. The first one is to construct a micro-local parametrix just 
as in the case of the heat equation, and then to use stationary phase method 
to evaluate the singular terms in the expansion. This appears, however, 
quite a cumbersome procedure even in the simple case of closed manifolds. 
Instead of going into detail, we give in the last section just an idea of several 
alternative methods. For more consistent exposition we refer to [73], [71] and 




(158) 



[75]. 
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7 Billiards 



Billiards are model systems that for many reasons have earned great popular- 
ity in the investigation of various problems in mathematical physics. Among 
their major advantages are their conceptual simplicity and the good behavior 
in the semiclassical limit, that makes them easy to quantize. Moreover, they 
are so far the best candidate in the attempt to quantizing chaos. 
There are many books and articles that for some reason incorporate this item 
- starting with [T3] , [32] and the appendix of pQ for the classical dynamics of 
a billiard, through the Berry articles [10] and [UJ (and many others) and ar- 
riving at some modern aspects, found in [3T], [77], [S] and [5B]. We specially 
refer to [72J for the elliptic case and to [53], [SZ] - for a thorough review of 
various mathematical aspects, concerning ergodicity as well. 

The notion of a (classical) billiard system is quite an intuitive one: given 
a connected domain M in M n with piecewise smooth boundary DM = UjiVj 
(iVj being smooth boundary arcs), the billiard dynamics is defined as the 
dynamics of a free particle, constrained in M and reflected according to the 



Hard ball map is constructed: given an initial point x £ M together with the 
initial direction (or equivalently a point £ in the unit ball bundle over M), 
we parameterize the unit geodesic 7^ passing from £ with some parameter 
r (usually chosen to be the distance, measured along the geodesic arc - we 
shall call it 'time' here), and then the 'transformed point' G>(£) we define as 
the position and normalized velocity of the same free particle after time r. 
We note that billiard trajectories can also be defined as the the trajectories 
along which the singularities of the wave equation are propagated. Most 
of them are simple n-bounce orbits (n > 2), governed by the Snell law. 
However if the trajectory intersects the boundary tangentially rather then 
transversally (forming a so-called 'gliding ray'), one has to deal with com- 
putational complications. Therefore, it is common practice to consider only 
convex domains for which all closed trajectories are n-bounce orbits plus 
closed geodesies on the boundary (it two dimensions there is only one such 
'glancing' trajectory that is the boundary itself). 

Another important class are the two-fold trajectories called also bouncing 
ball orbits. Each convex billiard has at least two of them, one being the 

75 according to which the angle of reflection is equal to the angle of incidence 




That is how the bil- 
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diameter. For this special type of orbits, the notion of ellipticity is reduced 
to the constraint that the length / 7 should satisfy either / 7 < min(p ,Pi) or 
max(po,Pi) < l-y < Po + pi, where po,i denote the radii of curvature at each 
of the two endpoints respectively. Alternatively, 7 is said to be hyperbolic if 
ly > po + pi or min(p ,Pi) < / 7 < max(p ,Pi)- In particular, if 7 is a local 
minimum diameter, then / 7 < p + pi and for local maxima we always have 

l-y > PO + Pi- 

In the elliptic case the eigenvalues of P 7 are of the form e ±ia , where the real 
parameter a 7^ kn is given by 



At the same time, the eigenvalues in the hyperbolic case, which are of the 
form e ±A , A G K are calculated with the help of the formula 



No matter the simple definition, billiards might become very complicated in 
terms of dynamics. Fortunately, there are special cases when they appear 
much easier to study, thanks to some exceptional features, such as integra- 
bility, ergodicity, or discrete symmetries. In this paper our main concern is 
the elliptic billiard table, which is a typical example of an integrable system 
(its integrability follows from the integrability of the geodesic flow of the 
ellipsoid, by a limiting procedure), but we also pay attention to the chaotic 
case, as it is rather conceptual. 

Polygonal billiards on the other hand are interesting as they illustrate the 
role of corners in the boundary. They can exhibit both integrable (triangle, 
rectangle) and chaotic motion. Nevertheless we omit this topic in order to 
be as concise as possible. 



7.1 Integrable and Chaotic Billiard Dynamics 

Now we investigate a little bit more thoroughly both types of billiards we are 
interested in. 

In the integrable case trajectories typically come in families, parameterized 
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Integrable systems 






i. h. n .i. !■■ i i Mi- 






Figure 4: Here we see plots of several typical billiard tables - the integrable circle, 
rectangle and ellipse, and the chaotic ones - Sinai, stadium and cardioid 

by KAM fori^l while in the latter they fill regions of the phase space every- 
where dense - the dynamics is said to be ergodic. 

We consider below the two most typical integrable billiard tables, shown in 
the picture - the rectangle and the ellipse (and its limiting case - the circle, 
usually studied separately). 

For the rectangular table integrability is obvious - after a suitable folding, the 
table is transformed to a torus and the billiard flow, up to a parametrization, 
coincides with the geodesic flow on the torus, which is integrable by definition. 
For the ellipse we may deduce the integrability from the integrability of the 
geodesic flow over the three-dimensional ellipsoid, which projects exactly on 
the billiard flow after shrinking one of the axes to zero size. This procedure 
can be shown to preserve the integrability. 

In order to show that the n— dimensional ellipsoid has integrable dynamics, 
76 the tori that foliate phase space, given by action-angle variables 
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Figure 5: The elliptic coordinates allow for the integration of the billiard table 
and the geodesic flow of the ellipsoid in the higher dimensional case. At the same 
time they describe most naturally the caustics - confocal ellipses and hyperbolae. 

one may introduce Jacobi elliptic coordinates A«, defined as the roots of the 
algebraic equation 



where are the semi-axes and Xk - the cartesian coordinates of x. However, 
by doing so one may attempt to show directly that the elliptic table is inte- 
grable. In the two-dimensional case these coordinates reduce to the sum and 
difference of the distances from a given point to the foci f\ and f 2 - So, let 
for example r(/i, f 2 ) = 2c, and ri(x) = r(x, fi) for each point x on the table. 
Then the elliptic coordinates are denned, according to the above definition 



In this way the whole table is knitted by an orthogonal net of confocal el- 
lipses £ = const and hyperbolae n = const, as shown on the picture. Now, 
using the fact that these are actually caustics for the closed billiard orbits and 
each orbits remains tangent to the same ellipse or hyperbola of the family, one 
may treat them as first integrals and thus show that the billiard is integrable. 

For chaotic billiards, in many cases the same limiting procedure can still be 
used. Planar convex tables with additional scatterers in the inside for exam- 




by 



£ = n + r 2 , n = n-r 2 
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ple0 may be regarded as projected hyperbolic surfaces. The latter are known 
to posses chaotic geodesic flow as seen from the Jacobi equation. Note also 
that the scatterer's boundary is concave from the inside of the table, which 
is reason enough for a chaotic motion. Other typical cases of chaotic tables 
are represented by the stadium billiard and the cardioid, shown above, as 
well as various Sinai billiards, that are usually constructed by removing a 
disc from the skin of the torus. In fact, a slightly modified Sinai billiard - 
a circular scatterer in a rectangular box, instead of a torus, would be the 
simplest known so far proof of the ergodicity of hard-sphere gasses - even for 
two particles, the system is chaotic, as the additional scatterer (representing 
the second particle) looks concave from the inside. 

In the case of chaotic billiards what one usually does in attempt to study 
the spectral statistics is to retrieve information about the mean level spac- 
ing, clustering etc., using ergodic theory. These estimates are often the only 
relevant data to predict and test by experiment. However in many cases they 
reach to amazingly rich results. 

It would be unfair not to mention the large class of billiards, known as 
pseudo-integrable. The pseudo-integrable systems in general are given by 
a hamiltonian 

H = H (p,q) + eH'(p,q) 

where H (p, q) is a completely integrable hamiltonian and e - a small param- 
eter. 

It is known from KAM theory, that for e close enough to zero, most tori 
survive and the set measure of totally destroyed tori grows with e, so the 
system is by no means ergodic, nor integrable for a generic H'(p, q). 



that arc somewhat natural generalization to the Sinai billiard 
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7.2 Billiard Quantization 



Roughly speaking, to quantize a billiard, means to study the wave (or Schroding 
dynamics instead of the geodesic (or billiard) flow. Due to the compactness 
of the domain in interest (the billiard table Q), the wave dynamics is reduced 
to the stationary equation 



However, since f2 is not anymore a closed manifold, one needs to impose suit- 
able boundary conditions on dQ. The most naturally appearing in quantum 
mechanical setting boundary conditions are homogenous Dirichlet. They are 
interpreted as infinite potential walls, surrounding the domain. Neumann 
boundary data is mostly applicable to problems concerning electrodynamics, 
acoustics and fluid dynamics. 

The Robin case is also relevant, as observed by Ballain and Bloch, in the 
investigation of a particular problem of nuclear physics. 

7.2.1 A Toy Model for the Toy Model 

The first explicit example to consider here will be the circular billiard table - 
one of the few that exhibit exact solutions. The dynamics is rather simple - it 
is governed by the conservation law p = const (all reflection angles are equal). 
Nevertheless, there is a conceptual difference between orbits with angles that 
are rational, respectively irrational sub-multiples of ir. The former close after 
n bounces while the latter fill a whole annuluJ^I. Note that even rational 
orbits are non-isolated, but come in families instead - we may rotate each 
rz-bounce orbit to an arbitrary angle E (0,2ir/n) that actually labels the 
family. In addition, all closed orbits are neutrally stable. 
The reason to consider this case in the first place is that the Helmholtz 
equation for the circle has an exact solution given by 



and the free Green function can be expressed in terms of modified Bessel 
functions as 



78 these are certainly projections of rational and irrational windings of the torus, denned 
by the integrability of the system - the inner radius of the torus projects onto the caustic 
for the corresponding family of orbits 




(159) 



ipim = ci m Ji(k m r) expiW, I E Z, m E Z/{0} 
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What one needs to do next is to adapt this kernel to the boundary conditions 
and then, using the classical Sohotski formula, express the spectral density 
in the large k limit. 

7.2.2 The Beauty of Ellipses 

Elliptic billiard tables posses some noticeable properties, in addition to inte- 
grability, that make them an interesting item in themselves. For example, it 
is a familiar result that a trajectory, passing trough one of the foci, passes 
trough the other after the first reflection. Moreover such trajectory asymp- 
totically tends to the major axis - the line, determined by the two foci. This 
is the stable bouncing-ball orbit. There is only one more two-fold trajectory 
- the minor axis, which is unstable. Note that these two lines are also the 
symmetry axes of the ellipse - it is preserved under reflection with respect 
to them. The fact that the table is parameterized by a caustic net is rather 
geometrical, even if we leave aside its crucial role for the integrability of the 
system. The two types of caustics determine two types of trajectories - the 
first, sometimes referred to as rotation-like trajectories are those, that never 
cross the segment between the foci - the ones that are tangent to confocal 
ellipses. Orbits in the second class are said to be of oscillatory type, as they 
bounce near the minor axis, each time crossing the segment between the foci 
and always touch the same confocal hyperbola. Note that according to our 
previous definitions these are exactly the elliptic and hyperbolic orbits in the 
billiard. We also make a point here that the fixed set of this billiard is clean, 
meaning that all fixed orbits are sub-manifolds, with their tangent spaces 
being the fixed set of the differential map, or, equivalently, the length func- 
tional is Bott-Morse on the free loop space. The whole chord space has the 
topology of Mobius strip, since closed orbits are periodic and invariant under 
time-reversal. When unfolded, this strip is seen to be divided into regions 
of elliptic and hyperbolic type of motion - note that due to the additional 
symmetries it looks rather specific. The separatrices are given by the two 
axes of the ellipse. 

As far as classical dynamics is concerned, we remind that the billiard map is 
integrable with the help of the elliptic coordinates, we previously introduced. 
If the boundary ellipse with foci, positioned at / 12 = (±e, 0), is defined by 
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Figure 6: The phase portrait of the elliptic billiard in s, p-coordinates: the two- 
bounce orbits may be regarded as degenerate cases of oscillatory and rotation-type 
motion, respectively. The stable one begins at p = 0, s = L/2, and the unstable 
(the minor axis) - at p = 0, s = 

the equation 

x 2 + zr—i = 1, < e < 1 
1 — t l 

then the elliptic coordinates n = coshp and £ = cos0 are introduced via the 
transformation 

(x, y) = (e cos cosh p, e sin sinh p) 

The confocal ellipses are the curves p = const, <fi e [0, 27r) and the hyperbolae 
- 4> — const, p E [0, arcosh(-)]. In these coordinates the classical hamiltonian 
of a free particle of unit mass is given by 

2e 2 (cosh 2 (p) -cos 2 ((f))) 

This hamiltonian is now most easily seen to be integrable - one just needs 
to multiply by the denominator and introduce the separation constant K to 
reduce the problem to a system of one-dimensional equations 

p 2 = 2Ee 2 cosh 2 p — K 

pi = K - 2Ee 2 cos 2 <p (160) 
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Then the two integrals of motion are the energy E and k 2 = K/2E. The 
above equations can be regarded as describing hamiltonian dynamical sys- 
tems with effective energy E p = —Ek 2 , E^ = Ek 2 and effective potentials 
V p = —Ee 2 cosh 2 p and = Ee 2 cos 2 4>, respectively. 

We note that the above representation allows also for exact quantization - 
the system is transformed to two identical Mathieu equations of the form 

(1 _ x i)f> _ x f + ( A _ cV)/ = 
with x = £ and x = r\ respectively. 

The eccentricity e of the ellipse is given in these coordinates by e = (cosh 2 p) _1 . 
If we consider the slope angle of the ellipse ip (the angle, between the coun- 
terclockwise oriented tangent and the at each point), it is obviously 
connected with the elliptic coordinates by the formula 

tan tjj = = — tanh p cos 4> 
ax 

and the radius of curvature is given correspondingly by 

ds e cosh p sinh p 



dip (cosh psin ijj + sinh p cos 2 i/j) 3 / 2 



In these coordinates one may find explicit expression for the caustics in the 
billiard as 

= p 2 - e 2 cos 2 jj(s) 
' 1 — e 2 cos 2 ip(s) 

In certain considerations, however, it is more convenient to work in different 
coordinates. Namely, let s denote the arc-length coordinate of the initial 
bounce of the orbit and p = cos 9 - the normalized momentum!^! Then the 
canonical symplectic form is given by the area element uo = sin 9 ds Ad8 and 
it is easy to see for example that the billiard map preserves it. 



7.3 Back to Semi-classics 

Various techniques have been involved so far in the semi-classical study of 
billiards. One can attempt to apply any of the methods discussed above as 

79 since the hamiltonian is homogenous function of p, it can always be re-scaled, restrict- 
ing the momentum to a unit ball sub-bundle of T*(dtt) 
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Figure 7: Standard plots of closed geodesies of different type for the elliptic bil- 
liard. For higher number of reflections caustics become quite visible. 



well as many others that have been invented for specific cases. 
For instance, one may go for the heat kernel expansion for the billiard ta- 
ble, considered as a manifold with boundary. Here are the first few heat 
kernel coefficients for a Dirichlet quantum billiard in terms of the boundary 
curvature k 



a = AreaiVt), a 1/2 = ~Jt; 



L(E) 



ai = ^ J Kda, a 3/2 = -gj J k 2 da 



4 f 3 j 37^ f 4 v/7T /" , 2 

° 2 = 315 y s K ' 05/2 = ~2i3~ y s K 2io y s ^ ^ 

where = j^da denotes the length of the boundary S = dQ. 

We note that for many purposes, as long as billiards are concerned, the heat 

trace approach is insufficient as it does not account for the dynamics in any 

explicit way - unlike the wave group, the heat semi-group is dissipative and 

could not be related directly to a conservative dynamical system. 

A more adequate approach in this direction would involve the wave trace, 

or semiclassical propagator traces. In the case of two-dimensional convex 

billiards for example, we have the two parts of the spectral density d(E) = 



107 



Figure 8: Plots of the probability density of different states for the elliptic table. 
The accumulation of intensity close to classical caustics is obvious. 



d(E) + d osc (E) in the form 

~ (^f^E^ (|v^ + 2,,-l/4) 

where U denotes the length of the i th orbit and Qj - the area of its band. The 
summation is performed over all primitive periodic orbits and their repeti- 
tions. 



7.3.1 Integrability vs. Ergodicity 

The problem of quantizing classically chaotic systems has drawn a lot of at- 
tention of both physicists and mathematicians. Unlike the integrable case, 
exact quantization here is inapplicable since uncertainty principle and dis- 
crete spectrum lead to somewhat regular behavior. The only way known so 
far for quantizing chaos is to consider instead the semi-classical limit h — > 



108 



of a given chaotic system. Only then the typically chaotic classical properties 
such as ergodicity or mixing are preserved 80 !. 

One of the fundamental systems in quantum mechanics - the helium atom 
(that is actually the three-body Coulomb problem) has been successfully 
quantized only with the help of the Gutzwiller trace formula, that is already 
a serious achievement. For detailed study on this problem we refer to [23] . 
Now let us turn our attention to the chaotic billiard tables, which are the 
main tool to model quantum chaos on the semiclassical level. Such systems 
differ from the integrable ones in several aspects. Firstly, their closed orbits 
are exponentially unstable and this usually leads to ergodicity in the phase 
portrait. Secondly, most orbits are isolated and pass trough almost every 
point of the 2N — 1 energy surface in phase space, due to the lack of first 
integrals. Nevertheless there are exceptional families of orbits, but these cer- 
tainly do not fill resonant A-tori as in the integrable case. 
There is significant difference between isolated and non-isolated orbits in 
terms of eigenvalue clustering as we shall see further. An adequate way to 
describe the first case is the Gutzwiller trace that has proved a fairly powerful 
tool. 

Nevertheless we follow here a slightly different approach. First we suggest an 
alternative description of quantum mechanics that is very popular, although 
not taught in the under-graduate course, and appears to be quite convenient 
for semi-classical considerations. Then we briefly compare the spectral statis- 
tics in the integrable and the chaotic case. 

The Wigner Function Approach The idea of the Wigner function for 
a quantum state if) is expanding in a basis of plain waves. More precisely, it 
is defined as the distribution 



The quantum-mechanical description of a system in terms of allows for 
representing in a convenient form basic physical quantities, such as the den- 
sity 



80 which means that the small h and large t limits do not commute, or the semi-classical 
estimate is not uniform in time 



W^(x,p) = (2nh)~ m [ 4){x - y/2)iP(x + y/2)e- lkM d n y 




109 



or in momentum representation 



\HP)\= J d m qW4q,p) 

but it also allows for a very 'nice-looking' formulation of the semiclassical 
description. Namely, we may insert the de Broigle wave approximation of ip 
in the definition of and compute the integral for small y. What we get 
is the form 

'd 2 s( q ,py 



w. 



det 



6(p-p(q,P)) 



dqidqj 

which, apart from other things, restricts the dynamics (in the roughest ap- 
proximation) to classical Lagrangian manifold, given by the Hamilton- Jacobi 
equation. 

In the integrable case we have a quantized torus, labeling functions by quan- 
tum numbers k — (k% . . . k m ) and 

W}(q,p) = ^S(I(p,q)-I k ) 

where the action variables are given by the WKB quantization condition 

h = h (k + ^) 

In the analogous situation concerning ergodic systems, we have similarly 

5(E-H(q,P)) 



J J d m pd m q5(E - H(q,p)) 



Using this approach one may find some qualitative differences between the 
two cases in terms of spectral statistics. We shall only enumerate some ba- 
sic results, concerning spectral clustering and degeneracies. A more detailed 
study on that topic may be found for example in [T3] . 

To begin with, we make a brief discussion on degeneracies. We shall be inter- 
ested not in degeneracies, coming from symmetries for the Hamiltonian, but 
only those, obtained by smooth varying of parameters in a family of 'simi- 
lar' Hamiltonians. Thorough investigation shows that for a 'generic' system, 
crossing of level curves demands varying at least two parameters (see [13] 
for details), so it is of codimension twd_j. This makes it highly unusual to 



81 in two dimensions these singularities are projections of the diabolic points at which 
cones in phase space meet 
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have such situation in ergodic or pseudo-integrable systems, although it is 
quite common for ergodic systems to have levels that become very close to 
each other and then repulse again (these may be viewed as hyperbolae on 
the touching cones, passing close to the singularity). 

However, in the case of torus-quantized systems, codimension one usually 
suffices (lattice planes always cross each other when perturbed), so degener- 
acy of such kind may be observed and a lot of examples are available. 
The level spacing S is another quantity that differs a lot between the ergodic 
and the integrable case. Let us consider the distribution P(S) defined as the 
differential probability for the spacing 



< AA > 



to lie in the interval (S, S + dS) and its asymptotic behavior as S — > 0. If it 
tends to zero, this can be thought of as effective 'repulsion' of levels, whereas 
if the limit of P is a non zero constant, then spectral levels become thicker 
and form a cluster at A. 

It has been proven that, broadly speaking, P ~ S k ^ 1 where k is the minimal 
codimension of the singular set of spectral degeneracies in the hamiltonian 
family, considered above. Thus, we have level repulsion for ergodic systems 
and clustering for those that exhibit torus quantization (integrable or pseudo- 
integrable). For integrable systems with more than one degree of freedom 
one has an even better estimate of the form P ~ e~ 5 . 

One more point to make here is that the mean level density in the Gutzwiller 
trace formula is of order h~ m so one might quite naturally expect to obtain 
it by smoothing d(E) over this range. However, this is not a good idea, since 
the oscillatory contribution has wavelength of order fi~ x which is much larger 
scale in dimensions higher than one. Therefore d and d osc should be treated 
separately. Let us focus for the time being on the latter. 
Note that for isolated orbits the amplitudes in d osc (E) oscillate with each rep- 
etition I for stable orbits and decay exponentially with / for unstable ones. 
In the integrable case, where orbits come in m — 1 parameter families, these 
amplitudes are shown to behave like l~^~ 

On the other hand each closed orbit contributes to the level density an os- 
cillation of wavelength given by 

\ j = T j (\)h/p 
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where Tj is the period and p - the momentum. According to this formula, 

closed orbits describe level clustering of order h, whereas the mean level 

spacing, as we already know, forms on the much finer scale h~ m . 

In the case of billiards, where the classical action along the j th orbit is given 

by 

Sj — — f%k£j 

we have an analogous rule for the wave-number range 

k = 2*K/tj 

Note also that in both cases each closed orbit does not give individual lev- 
els, but collective properties of the spectrunQ namely, clustering of range 
given by the above equations. Thus the determination of an individual level 
in the spectral density involves summation over infinitely many (or at least 
sufficiently many) closed orbits in order to generate a single delta function 
via interference of oscillatory terms. 

The sense of 'sufficiently many' varies from integrable to ergodic systems 
and for the former this number behaves like /j,-" 1 *™- 1 ) as J% _» o. It is even 
bigger in the ergodic case. For the Sinai billiard for example, as a typical 
representative, it has been evaluated to n ~ h~ 3 e a ^ h . 

In both cases there are much more practical methods for determining indi- 
vidual levels, but nevertheless the path summation is appreciable in certain 
cases - for instance, if we only want to know the 'smoothed' over a certain 
range AA level density, it suffices to include much less paths as the range 
increases. This can be of practical interest for acoustics, where certain reso- 
nances (in a concert hall for example) dominate the fine structure of spectra 
and thus introduce a natural 'smoothing' criterion. 

Now let us give some explicit examples of particular results obtained in the 
ergodic case. 

The first one is the Sinai billiard, that is a typical classically chaotic system. 
The smooth part of the spectral staircase function is given by 

M{E) = A(i _ ^) - ^[1 + V2/2 - r(2 - tt/4)] + § 

OZ7T 47T 9b 

82 this non-locality can be partly explained by the Fourier transform used in the deriva- 
tion of the Poisson type summation formulae and seems to be a major property of classical- 
to-quantum phenomena 
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where r is the radius of the scattering disk and the torus is taken to be 
double-folded unit square. The first term above stands for the area contri- 
bution, the second - the one of the perimeter, and the last represents the 
contribution from curvature and corners. 

A more thorough investigation would reveal the presence 'ghost' orbits - these 
are non-physical trajectories that pass through the obstacle (or, as in the case 
of non-convex tables, leave the interior of the billiard domain). One of the 
advantages of path summation in this case is that such orbits 'cancel' each 
other, after adding sufficiently many terms in the trace. For details we refer 
to [13]. 

For the case of Robin boundary conditions we refer to where the bound- 
ary function S is taken into account to obtain for the smooth part of the 
spectral density^ 



d(\) = Area{ ^^ 



27T 47T 



1 - 



+ 



whereas the oscillatory contribution comes with an additional phase (com- 
pared to the Dirichlet case) 

$ 7 ~ 2 arctan( ) 

where Pi = cos 9i are the tangential projections of the momenta in the points 
of reflection of the corresponding closed orbit 7. 



We leave this quite interesting topic of quantum chaology for the moment, 
realizing that there are many aspects and that we consider here only a few 
examples. For the most part we follow [13]. One may, however, decide to use 
more statistical methods, as it is done in [19J, or study more thoroughly the 
quantization of discrete maps, as shown in [23]. We refer to the latter for a 
thorough investigation on the matter. 

There is one more case we only mentioned in the previous section - the 
geodesic flow on a hyperbolic surface that also exhibits chaotic motion, but 
this time our research is armed with the more exact Selberg trace formula 
and the whole machinery of harmonic analysis on homogeneous spaces. For 

83 see the article for much more precise results 
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a review on this special case we refer to [17] and other articles by the same 
author. 
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8 Inverse Problems 



Solving an inverse spectral problem meas to determine geometrical data (met- 
rics, boundary conditions etc.) from spectral data. If isospectral deformation 
is impossible, then the metric or the boundary data is said to be spectrally 
rigid. If it is the unique geometry with this spectrum^!, then we call it spec- 
trally determined. 

There are only few positive inverse spectral results known so far and a vast 
number of counterexamples. For instance, it has been proved that the stan- 
dard metric on the sphere is spectrally determined only in dimensions up 
to six. The ellipses in the plane are known to be spectrally determined for 
the Robin boundary conditions, as we shall see in the following, but it is 
not known whether they are in the Dirichlet case. Flat tori are spectrally 
determined only in the low dimensions and in the higher-dimensional case 
this is so only locally - there are isospectral pairs that are not isometric, but 
no isospectral deformations are available, since each metric is determined in 
a sufficiently small neighborhood of the superspace. And in the end we know 
almost nothing about the hyperbolic manifolds. In dimension two we know 
that isospectral pairs contain only hyperbolic representatives, which cannot 
be said for sure in the generic case. 

However we have spectral determination within the corresponding class of 
simple real analytic surfaces of revolution with one critical distance from the 
axis and convex analytic plane domains with the symmetries of the ellipse 
and a bouncing ball orbit of fixed length. 

Spectral problems may be understood as revealing the relation between dy- 
namics, geometry and spectrum, hidden in the trace formulae. Note that in 
these formulae we cannot in general distinguish between different geodesies, 
but encounter just what is usually referred to as length spectrum. 
The length spectrum is defined as the set of all lengths of simple closed 
geodesies. Moreover, when multiplicities of geodesic lengths are also taken 
into account, we refer to this as the extended length spectrum. Since the 
length spectrum gives no information about which length belongs to which 
closed geodesic, in practice it is often preferable to deal with the so called 
marked length spectrum which classifies the closed geodesies in homotopy 

84 usually it is understood 'within the same class', meaning that each time we try to 
determine a metric or domain, we do it relative to other representatives of some restricted 
class of metrics and domains 
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types (in the case of closed manifolds) or in terms of rotational number^ 
(as long as convex domains are involved). 

It is curious to note that equality of the marked length spectra implies isom- 
etry at least for locally symmetric spaces and so does symplectic conjugacy 
between the geodesic flows. 



8.1 Basic Strategies 

How do we find a positive or a negative inverse spectral result? 
The negative one is certainly easier to prove, since one may just find a coun- 
terexample. For a positive result, however, a lot of work needs to be done. 
Perhaps the simplest inverse result is the spectral determination of Dirichlet 
discs. It uses only the two leading terms in the heat asymptotic expansion 
and the property that for any planar domain the inequality 

Area(n) < B^L 
holds and turns to equality if and only if Q is a disc. 

This case is, however, rather exceptional - there are many others in which the 
heat invariants remain constant under deformation which does not in general 
preserve the spectrum (for example the 'bumpy metrics' considered in [51] ) 
so one needs to resort to more subtle spectral data, such as the wave trace. 
To sum up, solving a generic inverse spectral problem is far from being trivial, 
but if there is any chance for it to be done, it would be best to have certain 
procedure to follow. The first step is usually defining as much spectral in- 
variants as possible (from heat or wave trace, zeta function etc). Next one 
needs to express these invariants in geometrical or dynamical terms, as in the 
case of the heat kernel coefficients. In the wave case we have the coefficients 
in the asymptotic expansion arising in an analogous way from singularities 
(or non- commutative residua). Under certain simplifying assumptions (for 
instance a simplicity of the length spectrum or analyticity of the domain) 
one may find a very 'nice' set of geometrical and dynamical spectral invari- 
ants, such as the length spectrum itself, or the Birkhoff normal form. But 
then comes the third crucial step when we need to completely determine the 

85 this quantity is denned as the ratio of the winding number and the number of reflec- 
tions 
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geometry (or dynamics) out of these invariants, and this is, first of all, rarely 
possible, and second - there is plenty of hard work to be done. 

Livsic Cohomology Let us now go a little bit more into detail. Up to 
first order in perturbation theory, we have for the variation of eigenvalues 

X k = (A0 fc ,0 fe ) 

where (pk are the unperturbed eigenvalues, and A - the variation of the 
laplacian, which is, generally speaking, a ^DO with a leading symbol g. 
The linearized problem suggests defining the space of variation of the metric 
(the tangent to the supers-pace of M at g) for which 

g = A fc = 

for all indices k and closed geodesies 7. At this stage a special class of 
co-homologies, the so-called Livsic cohomologies, are already involved. The 
Livsic cohomological problem studies whether each cocycle J f = V7 nec- 
essarily has a potential in the sense / = Xh(F) where is the generator of 
the geodesic hamiltonian. This is connected to inverse theory in the following 
way: let the variation of the metric tensor satisfy 

g = V7 

so that it also preserves the extended length spectrum. When the cohomol- 
ogy is trivial one has by all means the above representation of g and then 
one has to study the harmonic analysis of the underlying manifold in order 
to determine whether such representation is possible at all. For negatively 
curved surfaces the answer is 'no' and thus no iso-spectral deformations are 
available in this case. Therefore such manifolds are spectrally rigid. This 
technique works fairly well in the case of closed manifolds - especially when 
all geodesies are closed or in the case of integrable systems, foliated by flat 
tori - the Livsic equation then involves integration over these level sets and 
Fourier analysis suffices. However this approach is generally ill-defined or too 
complicated for billiards. 



.1 - 
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The Parametrix On the other hand one may choose to construct a parametrix 
of the wave equation in analogy with the heat kernel. However such construc- 
tion becomes in this case rather complicated^ even for simple geometries. 
Therefore we prefer to briefly review two other more up-to-date methods in 
the context of concrete problems they can be successfully applied to. 
For a more consistent exposition we refer to [71] and [75]. Among the few 
surveys on inverse problems, we adhere mostly to |73j and [5"T] . 



8.2 Quantum Birkhoff Normal Forms and Surfaces of 
Revolution 

Let us first note that the wave group is a quantization of the geodesic flow and 
certain correspondence-principle relations hold. For instance, if two lapla- 
cians are conjugated with respect to a unitary FIO, the geodesic flows are 
conjugated with respect to a homogeneous symplectic transformation. At 
least locally, isospectrality implies (for simple length spectrum) symplectic 
equivalence around closed orbits (equal Birkhoff normal forms). In the next 
paragraph we remind the notion of a classical normal forms, that is some- 
what central for the theory of dynamical systems, and then obtain formal 
quantization. 



Classical Birkhoff Normal Forms The so-called Birkhoff Normal Forms 
treat hamiltonians close to a stationary point (such as periodic orbits, that 
will be our main concern). Roughly speaking, putting an equation or a 
hamiltonian function into normal form means finding simple approximate 
expression (linear, polynomial etc.) 

Compared to the generic case, however, here there are certain principle diffi- 
culties due to the symplectic nature of the problem. First of all, the transfor- 
mation that brings H into normal form should be canonical, so usual cutting 
the Taylor series up to the quadratic term, as we are used to in the context of 
ODE close to equilibrium, is not sufficient. By doing so, in this case one may 
lose qualitative information - for instance the linearized equation may predict 
Lyapunov stability, that is observed in practice and can be easily destroyed if 
one takes a higher-order correction. If the system is close to integrability, it 

86 nevertheless it has been successfully used for obtaining certain inverse results before 
new schemes were developed 
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could be possible to find a set of action-angle variables for which its normal 
form is integrable up to the higher-order remainder. However its influence 
may accumulate over a sufficiently large time scale and may end up in a 
chaotic motion. 

There is one more major obstacle - the appearance of resonances. Suppose 
that one has managed to put H into normal form close to stable equilibrium. 
Then its quadratic term (for n degrees of freedom) will give an expression of 
the type 

1 n 

H~-J2co k r k (161) 

k=l 

where the canonical actions are expressed in terms of the new symplectic 
coordinates P and Q as tj = -^(Pj + Qf) and ujj are the corresponding fre- 
quencies. The system is said to be resonant of order s if there exists a 
linear combination of the frequencies with integer coefficients kj such that 
^2 \ kj\ = s. The main theorem, due to Birkhoff, states then, that if there 
are no resonances of order s and smaller, then H can be put in normal form 
which is polynomial of order s in the (P, Q)-variables (that is [s/2] in the 

r's) and the reminder is of order (9(t[~H). 



Quantization and relation to wave invariants The strategy of normal 
forms works for integrable dynamics and we illustrate it in the case of sim- 
ple real analytic surfaces of revolution. The first step is to use the Laplace 
spectrum to construct the quantum Birkhoff normal form of the wave group, 
that is essentially, expressing the operator y/A as a polynomial of the action 
operators I k and the tangential derivative along 7. Then, one hopes to re- 
constructs the metric from its normal form. After imposing some additional 
assumptions on the geodesic flow, in many cases this turns out to be possible, 
although the procedure is far non-trivial. Here the aim is only to give an idea 
of it, avoiding technicalities as much as possible. 

First of all one conjugates \fA into the maximal abelian algebra of fDO's, 
generated by the action variables around a closed non-degenerate geodesic 

and the tangential derivative D s := —i~r- Actually, this conjugation is per- 

ds 

formed on the normal bundle of 7 and the result is governed by the spectrum 
of P 1 . In the elliptic case for example, one has all the action variables in the 
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form of oscillators 

i k = a k (D 2 yk +y 2 k ) 

where a k are the exponent factors of the linear Poincare map spectrum and 
(s, Uk) are the Fermi normal coordinates around 7. 
In the real hyperbolic case one has alternatively 

I k = fi k (D yk y k + y k D yk ) 

Under some additional assumptions that we are going to specify later, it 
turns out that for a non-degenerate closed geodesic 7 there always exists a 
microlocal FIO, W : T*iV 7 /0 -> T^S 1 x R n ) which conjugates yfK in the 
form 

Wy^AW- 1 = D s + ^ + y Mlh__I?l3 (162) 

L z — ' D h s 

where p k are polynomials of degree k + 1 in the action variables and 2H a = 

Next, one hopes to establish a one-to-one correspondence between the normal 
form coefficients and the wave trace invariants. In order to do so, it is helpful 
to use residual calculus techniques to show that 

a 7)fe = Res z=0 TrD$ e itVA ' ^K~ z 

After putting ^/A into normal form one hopefully obtains the equality 

i a 

Vldet(I-P 7 )| 

where the $fc's form a special class of polynomials in the p/s and hence in 
the ij's, and the last factor to the right is a character of a metaplectic repre- 
sentation. The important thing we need to know about the above formula is 
that it is in our case invertible, which means that the Birkhoff normal form 
coefficients can be derived from the wave trace invariants. 



First Integrals The last step in the procedure consists in determining the 
metric form the normal form. This is actually deriving the metric from the 
hamiltonian function, which is quite clear conceptually, but demands some 
additional efforts. Therefore we restrict ourselves to the case of simple real 
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analytic surfaces of revolution with simple length spectra. 
The metric on such a surface may be written in the form 



ds 2 = dr 2 + a(rfdQ 2 



and impose the additional condition that there is precisely one critical point 
at which a" < that corresponds to equatorial geodesic. Moreover, we de- 
mand the Poincare map to be of 'twist' typer^l. 

Surfaces in this class have completely integrable geodesic flows and the corre- 
sponding hamiltonians posses global real analytic action-angle variables. The 
first action variable in our case obviously can be chosen to be the angular 
momentum I\ = pg and then the second takes the form 



with r± being the two extremals of r lying on the annulus determined by 
I\. From the first part of I2, one may retrieve that (see [H] for calculations) 
both Y2± \ a '( r ±( x )) \ an d J2± \ a '( r ±( x ))\ _1 are spectral invariant S and hence 
a'{r) is determined by the spectrum. But since we have a(0) = 0, the shape 
of the surface is spectrally determined. 

8.3 The Balain-Bloch Approach 

The idea of the method is essentially using exact identities for the Dirichlet 
or Neumann resolvent in terms of the free resolvent. It is far more explicit 
than the other two methods and rather convenient to apply for bounded 
domains. It has first been used by Balain and Bloch for the Robin problem 
in a three-dimensional ball. 

Suppose for the time being that we are dealing with a convex analytic region 
in the plane and consider an r-fold bouncing-ball orbit Y of length L 7 along 
the y-axis. Then denote by f± the Taylor coefficients of the boundary, 
considered as a function of one variable and calculated at both ends of the 

87 which means, roughly speaking, that its differential rotates all tangent vectors in the 
same direction 

88 note that the integral of interest is Abel transform and hence invertible, we have 
chosen x to denote a(r) 
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orbit. 

In this case we have the expression 

det(Id-P 7 ) = I%deth 

where h = Hess(/) and let h± denote the matrix entries of its inversd^I 
Then, under the assumption that the domain is symmetric with respect to 
reflections in the x-direction (so that we may drop the '±' subscript) and 
7, 7 _1 are the only geodesies whose length is L 7 , one may express the wave 
invariants at 7 as polynomials of the /^'s. 

In order to see why this is so, let us first consider the boundary-value problem 
for the wave equation and construct a regularized resolvent R p using a cut-off 
p around the geodesic length L 7 . Under the above assumptions, this resolvent 
admits a complete asymptotic expansion of the form 

00 

Tr I n R p {k + it) ~ e {ik ~ T)L ~< ^ (a yj + a^k^, k -> 00 (163) 

i=i 

which is essentially the same as the wave trace asymptotic at t = L 7 . In order 
to find it, one reduces the problem to the boundary, using exact identities 
from potential theory. The key formula in the so-called Fredholm- Neumann 
reduction yields 

R n (k+ir) = R (k+iT)-Di(k+iT)o(Id+N(k+iT))- 1 oBoS£ t (k+iT) (164) 

where Ro is the free-space resolvent in R 2 and B : H S (Q) -> H"- 1 l 2 {dO) is the 
restriction to the boundary. The single and double-layer potential operators 
SI, D£ : H s (dtt) -> H*+ 1/2 (Q) are constructed as in (JH6J, with the help 
of the free-space Green function Go and its normal derivative respectively. 
Note that here Si and D£ are considered as Fourier integral operators in this 
context (S£ t denoting the transpose), and finally N(k + it) stands for the 
the Fredholm resolvent, that is the boundary integral operator, induced by 
D£ G \I /_1 (9fi)(see [73J). It is classical that Id+iV(A;+ZT) is an automorphism 
of the Sobolev space H s (dQ). 

The above formula is in close analogy to the Grushin result, that relates the 
trace of the wave group to traces, concerning the monodromy operator (see 
[73] for example). 

89 note that one has in addition h l £ — h l _T 
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Following this analogy a bit farther, under the assumption that L 7 is the 
only length in supp/5 one obtains 

/ P(k- A) — lndet(Id + N(k + iT))d\~\2a i:j k- j (165) 
Jr d\ 

where the the wave invariants at 7 and det denotes the usual Fred- 

holm determinant. 

Next one examines semi-classically N(k + ir) and the geometric series expan- 
sion for (Id + N(k + ir))^ 1 . The former clearly has a WKB approximation 
away from the diagonal, where it has singularity of order —1, whichl^l rep- 
resents it by a FIO with phase equal to the boundary distance function 91 ! 
d((f>,4>') = \q(4>) — q{4>') |. Then substituting the explicit form of the free 
Green function in the plane into the definition of N, one obtains for its 
integral kernel 

N(k + iT,q((j)),q((j)')) = -2(A; + ^r)iJ^ 1) (A; + ^rrf(0,0 , ))cosa(0,0 , ) (166) 



where a(<f>,(j>') = Z(q((j)) — q((j)'),u q ^\) with being the unit normal at 
the boundary point q(<f)). Note that this is exactly the angle of incidence 
at q(4>) for a billiard trajectory passing through q((f>) and q(4>'). Then keep- 
ing in mind that for convex domains the boundary distance function d(4>, <$') 
generates the billiard dynamics, it is natural to regard N(k + it) as a global 
quantization of the billiard ball map. 

Due to the singularity at the diagonal it is not possible to apply directly the 
stationary-phase approximation scheme and we need some kind of regular- 
ization. So the next step is to perform (at least formally) a finite geometric 
series operator expansion with a remainder. After introducing a suitable 
cut-off, one separates N(k + it) into a tangential and transversal part. The 
former is irrelevant to the semi-classical quantization of the billiard map and 
behaves like an Airy operator, classically associated to gliding (or 'creeping') 
rays at the boundary. Its contribution turns out to be possible to factor out, 
which is one of the main features of the Balain-Bloch approach. Another 
feature is the possibility to introduce a cut-off at 7 that ensures the good 



90 it is actually order —2 in the two-dimensional case 

91 from the off-diagonal asymptotic of Go it follows that as long as \q(4>) — ?(</>') I — l^| 1_e 
for some e < 1, N(k + it) is a semi-classical FIO with a phase function given by the 
boundary distance function 
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behavior of the estimate. 

After having done all this, one ends up with oscillatory integral series which 
demand a stationary phase method to evaluate. Therefore a kind of Feyn- 
man graph technique 92 ! is incorporated to determine all the coefficients. 
Calculations show that the coefficients f±^(0) and /± (0) of the Taylor 
expansion for the boundar}@ appear in the stationary-phase expansion for 
the first time in the term of order k~i +l . In the presence of a symmetry axis, 
the leading part of this term has the fornF*! 

a 7J _ 1 ~2rL 7 (/ i 11 )V (2i) (0) + ... 

where the dots stand for lower order derivatives. And in addition, if the 
domain has the symmetries of the ellipse, since the odd-dimensional coef- 
ficients disappear, one may iteratively determine the /' 2 ^'s and hence the 
whole geometry. This gives a strong result for analytic convex plane do- 
mains with two symmetries - they can be determined by the wave invariants 
at a bouncing-ball orbit along the second symmetric axis. 

8.4 Variation on a Theme by Guillemin and Melrose 

Actually, although we revise it in the context of the Balain-Bloch strategy, 
the results obtained in [36] and [37], rely on the construction of a microlocal 
paramatrix. 

The main result in [36] is that the spectrum of the Robin problem in the 
interior of the ellipse, together with the pure Neumann spectrum completely 
determine the function S, taking part in the boundary condition. More 
precisely, the difference of the wave traces for the modified and pure Neumann 
case respectively: 

X = Tr cos t-sJ A s — Tr cos r y/Ao 

is shown to be tyDO of order —1 (the most singular parts cancel each other) 
which admits an asymptotic expansion in fractional powers of t — T at t = T, 
around a closed geodesic with period T and leading symbol 

M*)~£ff (167) 

j J 

92 a diagrammatic technique used as a means of computing integrals appearing as quan- 
tum corrections in perturbation theory 

93 considered as a graph of a function of one variable 
94 recall the notations from the beginning of the paragraph 
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where summation is taken over all points of reflection at the boundary. 
Guillemin and Melrose menaged to rewrite this spectral invariant as an in- 
tegral over the set of closed orbit 9^1 with period T ( or below denotes the 
induced invariant measure) 

Mx)= / ~^da T (168) 
Jy t sin 8 

Using some simple arguments from the geometry of the ellipse, the authors 
rewrite the above in the form 

2?r S(s)h(s)ds 
^2M{s) + Z-b 

where the function Z actually defines the confocal caustics by the equation 

x 2 y 2 



e + Z Z 

(e = a — b as before) and M(s) denotes the function 

, , a + b ab 

M(s) = : 

2 a + b — e cos 2s 

Moreover we have 

and the rest follows from considerations concerning the invariance of the 
integral measure. More precisely, the volume form u = sm8d8 A h(s)ds can 
be written also as 

uj = dZ A f(s)ds 

where the second multiplier to the right is the invariant Leray form on the 
curves Z = const, that is exactly the one we need. We have from the above 

98 

f(s) = sin8—h(s) 

and easily find 

BZ 

— = -2y/(Z-b)(2M-Z + b) 



95 which is a smoothly imbedded manifold since the billiard map is clean and its action 
is conjugate to rotation 
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Using the symmetry of the ellipse D = Z 2 <S> Z 2 , one may write the above 
integral as a sum of four components and the function S is determined at the 
endpoints of each interval, that could be easily checked, expanding over two 
'special' orbits - the major and minor axes. Moreover, since the remaining 
factor in the integral is strictly monotonous in each of the intervals, one sees 
that the boundary function is determined itself. 

Evaluating the significance of the result of Guillemin and Melrose, we should 
admit that it is not very likely that in practice anyone will ever desperately 
need to 'guess' the elasticity of the boundary of a vibrating elliptical region 
by its spectrum, but this is one of the very few exactly solved inverse results 
that makes it of particular interest for further investigation. The articles [55] 
and [77] provide some insight about this. 

In this section we offer an easy way to obtain explicit formulae for the spectral 
variations with respect to varying the boundary function. It is worth saying 
that the method is completely identical to the one offered in the second 
chapter of [52] , where it was applied to the pure Neumann case only. As 
a matter of fact the same investigation has been cited in [77] , but using a 
different method, and obtaining, as it seems, a different result, which we 
consider a minor misunderstanding. 
To start with, we regard the spectral problem 

(A + A)0(x)=O xett (170) 

in a plane domain Q , enclosed by the ellipse X = dQ. 

Let 0o be a solution satisfying the homogeneous Robin boundary condition 

^- + S (s))^o(s) = S6S (171) 



dn 

and let Go denote the corresponding Green function. By definition we have 

(A + A) G° x (x, y) = 5{x-y) (172) 
(jL + S (x)^G° x (x,y) = 0, xeS 

Now let us perturb the boundary condition with a small correction 

S ^S = S + eS' + 0(e 2 ) 
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and denote with and G the solution and the Green function for S respec- 
tively. 

Then we take the spectral equation for G° x and multiply it by 0, and the one 
for we multiply by After subtracting the latter from the former and 
integrating over the Q we get 

0(x) = f [<f>(y)AG x (x,y)-G o x (x,y)A(f>(y)]du J (173) 

Applying the Green theorem to the above leads to 

4>(x) =e J S'G° x <f>dcr (174) 

from which is determined. 

In order to find A' we choose to work in the basis spanned by the eigenfunc- 
tions 0o of the non-perturbed system. We expand in this basis 

= ^c Q Q o (175) 

and take the usual spectral expansion of the Green function 

e!M = E^ ("<» 

a k 

which, substituted in (I174p . together with (I175p . leads to 

E^ = e E^;T~A^ V/i (177) 

where S' are the L^E) matrix entries of the boundary function derivative 
in the sense that 

S '^= I S'(o)(p^o)(p u (a)do 
Jan 

Then (11771) turns into a system of homogenous matric equations (for all //) 
for the unknown coefficients c v . More precisely 

(A-Aj) C/i = e XX^ ( 178 ) 

Now we restrict ourselves to the solution M of the perturbed problem, cor- 
responding to eigenvalue A M . If we would like to solve the integral equation 
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(11741) with the method of consequent approximations, we need to substitute 
as a first 'guess' <pQ under the integral sign. Then, up to first order we have 
c v = 0, v ^ fi and certainly — Then, this first approximation gives 

A,-Aj = eS; M + 0(e 2 ) 

meaning that the spectrum is being shifted up to first order in e namely by 
the diagonal element of the boundary function correction, exactly as in the 
standard stationary perturbation theory, only with the difference that this 
time the role of 'potential' is played by a single layer distribution instead of 
a volume term. 

In more explicit terms, we have for the first correction to the Robin spectrum 
the expression 

K = f S'(a)\<p»\ 2 da (179) 

The analogy to the old quantum mechanical perturbation-theory problem 
spreads even farther and allows for writing the higher-order terms as well. 
For instance one has 

A ~ e To 

and so forth. 

From the above considerations it also becomes clear how to vary the trace of 
the wave propagator 

TH7(r) = Trexp (iry/E) = ^V^^ 



According to the Guillemin and Melrose's result, we hav 

8a L \ T=T = e [ 4^rda T (180) 
Jy T sm6 l (s) 

According to [36] the dependence of the wave trace on the boundary function 
appears only in the second term, so it is not sensitive to varying S until the 
third term. Or, in other words, the ^>DO corresponding to the variation, is 
of order —2. 



96 note that the classical momenta p — cos 9 and thus the denominator sin 6* = yl — p 2 
are not altered by changing the boundary condition, nor is the integral measure 
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We note once more that the results obtained by Zelditch in [77] (see page 14) 
are slightly different - for instance he obtains a variation of the wave trace of 



the formal 



STtU(t) =ir22 e 



and this is used (again in cooperation with the result of Guillemin and Mel- 
rose) as a mens of evaluation of boundary trace^^ such as 

\ k <\ 

for the n-dimensional case. Note that this power dependence suffers slight 
modification according to the present considerations. 

Now let us see what happens if we formally differentiate directly the trace 
( 11801) . According to the simple rule for composite functions we have 

4 TiU(t) = it V e irV ^ (181) 

de\e=o V ; ^2^Y V 

Substitution into the earlier obtained form of \ v yields 



T Trf/(r) = - / ^2^S'(s)\r (s)\ 2 da (182) 

Note also that from the above expression one may easily find the first terms 
in the power series expansion for the wave trace perturbation. Namely 

|,„ „ T *M = ixC(l/2) - £(-1) - % E A„VX + £(-2) + . . . 

Then using again the representation of the zeta function as a Mellin trans- 
form of the heat kernel, we obtain after cancelation of individual terms, an 
expression for the linear term 



d 

de\e=o 



TrU(r) ~ J™ dt (2S't- 1 ' 2 + -^S' + 2SS') + 0(Vt)^j , 



T 



97 that obviously contradicts with our result 

98 the definition of the boundary trace of a function in the context of a certain boundary- 
value problem is f b = (I — B)f\gn, where B is the corresponding to the condition boundary 
operator 
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and may also fight for expressing higher-order terms in a similar manner. 
To the contrary, one may express ^-invariants in terms of derivatives of wave 
trace coefficients. 

Note also that the above formulae hold in the strict sense only away from the 
singularities, present due to periodic orbits, that is for r/T 7 . However, one 
would like to use it together with (I180p in order to obtain a formal expression 
for the first derivative of a^. 
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